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ABSTRACT 

Statistics  for  the  Global  Applicability  of  Retrieval  Algorithms 
for  Geophysical  Parameters  Using  the  SSM/I.  (August  1995) 
John  Francis  Polander,  B.S.,  Texas  A&M  University 
Chair  of  Advisory  Committee:  Dr.  Thomas  T.  Wilheit 


Brightness  temperatures  (Tb),  measured  by  the  SSM/I  (Special  Sensor 
Microwave  Imager)  at  85.5  GHz  (vertical  and  horizontal  polarization)  during 
August  -  November  1987  and  January  -  March  1993,  are  examined  over  the  Earth’s 
oceans  between  60°  N  and  60°  S.  Collected  data  are  placed  into  regions  of 
dimension  5°  latitude  by  5°  longitude.  For  each  region,  an  absorption  and  scattering 
regime  is  identified  based  upon  the  frequency  of  occurrence  between  the  average 
brightness  temperature  ((Tegsv  +  Tb85h)/2)  and  the  polarization  difference  (Tegsv  - 
Tb85h)  a  global  distribution  of  sensitivities  is  then  constructed  which  allows  us  to 
determine  how  often  our  new  and  improved  algorithms  can  be  used  for  computing 
various  geophysical  parameters;  cloud  liquid  water  (CLW),  precipitable  water 
(PW),  sea  surface  temperature  (SST),  and  sea  surface  wind  speed  (WS).  Results 
indicate  the  algorithms  can  be  applied  75-90  %  of  the  time  over  the  majority  of  the 
Earth’s  oceans.  These  algorithms  will  serve  as  a  simulation  study  for  the 
performance  of  the  Multi-frequency  Imaging  Microwave  Radiometer  (MIMR) 
which  is  scheduled  for  launch  in  the  year  2000. 

Additionally,  a  quantitative  assessment  of  CLW  distributions  for 
nonprecipitating  clouds  is  presented  over  selected  oceanic  regions.  This  gives 
insight  into  the  amount  and  variability  of  CLW  present  in  clouds  within  the 
radiometer’s  field  of  view.  The  results  indicate  little  variability  in  subtropical 
regions  with  a  median  value  around  .05  kg/m^.  Increased  variability  and  higher 
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median  values  (.07-.  1 1  kg/  m  )  are  characteristic  of  the  remaining  oceanic  regions. 
The  only  exceptions  are  along  the  ITCZ,  SPCZ,  and  monsoon  regions,  where 
frequent  precipitation  and  dense  clouds  render  the  atmosphere  sufficiently  opaque, 
leading  to  unreliable  results. 
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CHAPTER  I 

INTRODUCTION  AND  BACKGROUND 

1.  Introduction 

Water  vapor  and  liquid  water  are  strong  absorbers  and  emitters  of  radiation.  Ice, 
on  the  other  hand,  can  be  a  strong  scatterer  of  electromagnetic  energy  in  the  visible 
and  microwave  wavelengths.  Ice  also  absorbs  and  emits  radiation  in  the  infrared 
spectrum.  Gmnn  and  East  (1954)  demonstrated  the  dominance  of  scattering  versus 
absorption  of  radiation  by  ice  particles  at  a  given  wavelength.  The  theory  works  well 
and  is  relatively  simple  for  Rayleigh  scatterers  when  rj^e  «  K  where  rjcg  is  the 
radius  of  ice  particles  and  X  is  the  wavelength.  Such  is  the  case  for  85.5  GHz.  It  is 
also  applicable  when  rjce  >  X,  which  includes  the  Mie  scattering  regime  and  geometric 
optics.  In  those  instances,  the  mathematical  computations  become  somewhat  more 
cumbersome,  but  the  principal  results  are  unchanged. 

With  the  advent  of  computer  and  space  age  technology,  research  scientists  are 
now  capable  of  retrieving  vast  quantities  of  data  from  satellites  carrying  radiometric 
instruments  such  as  the  Special  Sensor  Microwave/Imager  (SSM/I).  These 
instruments  utilize  passive  remote  sensing  as  a  means  of  measuring  the  upwelling 
radiation  from  the  Earth’s  surface  and  atmosphere  and  relating  it  in  terms  of 
brightness  temperatures  (Tb).  In  particular,  microwave  radiometry  provides 
information  in  almost  all  types  of  weather  regimes  as  opposed  to  infrared  and  visible 
(Filiberti  et  al.  1994).  Because  microwave  radiation  is  absorbed  10  to  100  times 
more  by  liquid  water  clouds  than  by  ice  clouds,  we  can  deduce  cloud  liquid  water 
content  assuming  the  effects  of  water  vapor  are  taken  into  account  (Alishouse  et  al. 
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1990). 

The  primary  focus  of  this  study  is  to  determine  how  often  and  over  how  much  of 
the  Earth’s  ocean  surface  we  will  be  able  to  accurately  measure  various  geophysical 
parameters  using  the  Multi-frequency  Imaging  Microwave  Radiometer  (MIMR). 
These  parameters  include  sea  surface  temperature  (SST),  surface  wind  speed  (WS), 
cloud  liquid  water  (CLW),  and  precipitable  water  (PW)  sometimes  referred  to  as 
columnar  water  vapor  or  total  integrated  water  vapor. 

Obtaining  accurate  retrievals  of  SST,  WS,  CLW,  and  PW  is  extremely  important 
since  conventional  observations  over  the  oceans  are  very  sparse  and  in  most  places 
nonexistent.  Passive  microwave  radiometry  provides  the  means  to  access  data  from 
otherwise  inaccessible  oceanic  regions. 

The  major  drawback  to  obtaining  accurate  retrievals  occurs  in  the  presence  of 
precipitation.  Regions  containing  significant  precipitation  or  dense  clouds  render  the 
retrievals  unreliable  (Goodberlet  et  al.  1989).  Conversely,  in  the  presence  of  thin 
clouds  or  clear  skies,  the  retrievals  can  be  made  to  produce  better  results  (Nelkin 
1992).  Therefore,  the  elimination  of  data  associated  with  rain  and  knowing  where  the 
clouds  are  thin  are  of  prime  importance  and  somewhat  problematic. 

To  solve  this  problem,  a  filter  is  created  based  on  85-90  GHz  polarization 
differences  (Tegsv  -  Tbssh)  to  eliminate  observations  which  likely  contain  rain  and,  if 
desired,  excessively  thick  clouds.  The  remaining  observations  are  associated  with  a 
range  of  atmospheric  opacities  which  are  described  in  terms  of  their  85  GHz 
polarization  differences.  A  large  ATg  (typically  >  25  K)  indicates  a  relatively  clear 
atmosphere  (thin  or  no  clouds)  while  a  small  ATb  (<  10  K)  is  associated  with  a  nearly 
opaque  atmosphere  (thick  clouds  and  possible  rain). 

The  filter  is  also  used  to  select  algorithms  and  to  determine  how  often  the  more 
sensitive  algorithms,  for  a  given  parameter,  can  be  used.  For  example,  selecting  an 
algorithm  to  retrieve  the  value  of  a  specific  parameter  (i.e.  SST)  involves  several  steps. 
First,  all  the  85-90  GHz  observations,  within  the  radiometer’s  lowest  frequency 
channel  (6.8  GHz  in  this  example)  field  of  view  (FOV),  are  collected.  Then,  the  pair 
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(Tb85v  .  Tb85h)  of  observations  with  the  smallest  ATb  in  that  FOV  is  identified.  This 
pair  is  associated  with  the  most  problematic  part  of  the  FOV:  the  region  with  the  most 
opaque  atmosphere. 

After  identifying  the  pair  of  observations  with  the  smallest  ATb  >  an  algorithm  can 
be  selected  that  is  consistent  with  this  ATb  value.  Specifically,  this  will  be  the  most 
sensitive  algorithm  than  can  be  safely  used  while  remaining  compatible  with  the  most 
opaque  portion  of  the  atmosphere  within  the  FOV. 

Finally,  some  understanding  can  be  gained  regarding  the  ATb  distributions  in 
terms  of  the  amount  and  distribution  of  cloud  liquid  water  in  nonprecipitating  clouds. 

2.  SSM/I  and  MEMR  descriptions 

The  first  launch  of  the  SSM/I  (July  1987)  was  aboard  the  Defense  Meteorological 
Satellite  Program  (DMSP)  satellite  designated  F-8.  Three  additional  instruments  are 
currently  in  operation  on  DMSP  satellites  designated  F-10,  F-11,  and  F-13.  A  fifth 
SSM/I  aboard  the  recently  launched  DMSP  F-12  platform  will  soon  provide 
additional  data.  However,  technical  difficulties  are  currently  precluding  it  from 
operational  use. 

The  SSM/I’s  operate  in  a  near-polar,  sun-synchronous,  circular  orbit  at  an  altitude 
of  833  km.  Their  orbital  inclination  is  98.8°  (i.e.,  81.2°  retrograde)  with  an  orbital 
period  of  102  minutes.  The  conical  scan  from  the  SSM/I  covers  a  swath  of 
approximately  1400  km.  Depending  on  the  frequency,  the  instantaneous  spatial 
resolution  varies  from  69  x  43  km  at  19.35  GHz  to  15  x  13  km  at  85.5  GHz. 

The  SSM/I  measures  the  upwelling  radiation  from  the  Earth’s  surface  and 
atmosphere.  If  scattering  is  not  considered  and  horizontal  homogeneity  is  assumed 
then  by  integrating  the  Equation  of  Radiative  Transfer  (ERT)  and  using  the  Rayleigh- 
Jeans  approximation,  a  brightness  temperature  (Tb)  can  be  expressed  as: 
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}  dh' 

T(a,b)  -  Jy(h')--^ 

a  COS(e) 

Tb  is  the  brightness  temperature  measured  by  the  SSM/I  and  exiting  the 
atmosphere  at  the  incidence  angle  0  (average  of  53.1°  for  the  SSM/I).  T(h)  and  y  (h) 
are  the  thermometric  temperature  and  absorption  coefficient  at  height  h.  The 
reflectivity  and  earth’s  surface  temperature  are  r  and  respectively.  Applying  a 
derivation  of  Kirchoff  s  Law  (reflectivity  (r^)  +  emissivity  (e^)  =  1),  the  quantity  (1-r) 
denotes  the  surface  emissivity  while  T^b  is  the  brightness  temperature  of  the  cosmic 
background  (2.7  K).  The  absorption  coefficient  (y)  times  the  layer  thickness  is  the 
optical  depth  (x)  between  points  a  and  b. 

Term  A  denotes  the  upwelling  radiation  from  each  incremental  atmospheric  layer 
and  attenuated  by  intervening  layers.  Term  B  denotes  the  surface  emission  attenuated 
by  the  atmosphere.  Terms  C  and  D  result  from  the  radiation  reflected  off  the  surface. 
Term  C  denotes  the  downwelling  emission  attenuated  by  the  segment  of  the 
atmosphere  between  the  emitting  layer  and  the  surface  and  by  the  entire  atmosphere 
on  the  way  back  up  to  the  sensor.  Term  D  denotes  the  attenuation  of  the  cosmic 
background  (2.7  K)  from  its  downward  and  upward  pass  through  the  atmosphere  and 
reflection  off  the  surface  (Wilheit  et  al.  1994). 

The  SSM/I  is  a  passive  microwave  radiometer  which  measures  brightness 
temperatures  at  four  frequencies:  19.35,  22.235,  37.0,  and  85.5  GHz  in  both  horizontal 
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(H)  and  vertical  (V)  polarization  (from  this  point  on,  frequencies  and  their 
polarizations  will  be  referred  to  by  their  integer  value  followed  by  a  “V”  or  “H”  if 
applicable).  The  exception  is  at  22  GHz  where  only  the  vertical  polarization  is 
measured.  This  totals  seven  separate  channels. 

The  SSM/I  spins  about  the  vertical  axis  with  a  rotation  period  of  1.9  seconds. 
Observations  are  taken  in  a  102.4°  portion  of  the  rotation  with  the  SSM/I  pointed  in 
the  aft  direction  aboard  the  F8  satellite  (forward  direction  for  all  others).  The  102.4° 
arc  is  centered  on  the  satellite  subtrack  and  covers  a  1394  km  wide  swath  given  an 
orbital  height  of  833  km.  During  each  scan,  within  the  102.4°  arc,  both  85  GHz 
channels  are  sampled  128  times.  The  distance  between  samples  along  and  between 
scans  is  approximately  12.5  km.  The  integration  period  for  a  single  sample  is  3.89 
ms.  This  results  in  the  85  GHz  channels  having,  at  the  Earth’s  surface,  an  effective 
resolution  of  approximately  1 5  km. 

The  channels  at  the  lower  three  frequencies  (19,  22,  37)  are  sampled  every  other 
scan.  The  scans  where  all  channels  are  sampled  are  designated  ‘A-scans’.  The  scans 
which  contain  only  the  85  GHz  observations  are  designated  ‘B-scans’.  The  lower 
channels  are  sampled  64  times  along  an  A-scan  (vs.  128  for  85  GHz)  with  an 
integration  period  of  7.95  ms. 

Channel  resolution  varies  as  a  fimction  of  frequency,  antenna  size,  and  altitude. 
Taking  into  account  the  integration  period,  the  effective  resolutions  of  the  lower  three 
frequencies  are  approximately  69  km  (along-track)  x  50  km  (cross-track)  at  19  GHz, 
60  km  X  47  km  at  22  GHz,  and  37  km  x  38  km  at  37  GHz. 

Additional  details  on  the  sensor  design  of  the  SSM/I  and  specifics  for  computing 
brightness  temperatures  are  located  in  the  User 's  Marmal-SSM/I  Antenna  Temperature 
Tapes  (Wentz  1988)  and  the  SSM/I  User 's  Guide  (Hollinger  et  al.  1987).  Appendix  A 
is  a  schematic  showing  how  the  SSM/I  views  the  polarizations. 

The  Multi-frequency  Imaging  Microwave  Radiometer  (MIMR)  is  a  planned 
future  passive  microwave  instrument  that  will  measure  brightness  temperatures  over 
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TABLE  1.  Intercomparison  of  MIMR  and  SSM/I  performance.  (Values 


taken  from  the  MIMR  interim  report,  15  November  1994,  draft  0). 


Performances 

MIMR 

MIMR 

SSM/I 

EOS-PM 

METOP 

Altitude  (kms) 

705 

824 

833 

Incidence  angle 

50 

55 

53.1 

Channels  (GHz) 

6.8/10.65/18.7 

Same  as 

19.35/ 

23.8/36.5/89.0 

EOS-PM 

22.235 

37.0/85.5 

Sensitivity  (K) 

@  6.8  GHz 

0.2 

Same  as 

- 

@  19  GHz 

0.5 

EOS-PM 

0.8 

@89  GHz 

0.7 

1.1 

Footprint  size 

(kms) 

@  6.8  GHz 

60 

85 

- 

@  19  GHz 

22 

31 

69 

@89  GHz 

4.9 

6.8 

15 

Swath  (kms) 

1230 

1394 

Antenna  size  (m) 

1.6 

1.6 

0.61 
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12  channels  (6  frequencies)  ranging  from  6.8  to  89  GHz.  It  is  being  considered  for 
launch  at  the  turn  of  the  century  on  two  separate  platforms;  the  METOP  and  US  EOS- 
PM  satellites.  See  Table  1  for  a  comparison  of  SSM/I  and  MIMR.  performance 
features. 

3.  Previous  work 

Recent  research  has  been  conducted  for  developing  algorithms  to  estimate  sea 
surface  and  atmospheric  parameters,  specifically  SST,  WS,  CLW,  and  PW.  Several 
studies  were  conducted  which  discussed  the  development  and  the  improvement  of 
algorithms  used  to  calculate  WS  (Wilheit  1979;  Goodberlet  et  al.  1989)  and  both  WS 
and  direction  simultaneously  (Wentz  1992).  Microwave  emission  from  the 
ocean  varies  according  to  the  sea  surface  state,  i.e.,  its  roughness  and  foam 
coverage.  Calm  seas  are  characterized  by  a  relatively  smooth  surface  which  is  highly 
reflective  (low  emissivity)  and  is  associated  with  a  highly  polarized  emission, 
assuming  the  satellite’s  view  angle  is  far  away  from  nadir  (as  with  the  SSM/I).  The 
Fresnel  relations  from  electromagnetic  theory  (Griffiths  1989)  become  applicable  and 
state  that  the  reflectivity  for  horizontally  polarized  radiation  is  greater  than  for  the 
vertical;  therefore,  the  horizontally  polarized  brightness  temperature  is  always  less 
than  the  vertical.  If  the  sea  surface  becomes  rough,  the  polarizations  become  mixed 
so  that  the  difference  (ATg )  decreases  (Wentz  1992). 

Sea  foam  and  roughness  are  mechanisms  that  contribute  to  the  variability  of 
surface  emissivity  and  reflectivity.  Wilheit  (1979)  used  a  quantitative  model  to 
describe  their  effects  on  the  microwave  emissivity  of  the  ocean  surface.  He  compared 
his  model  with  observations  obtained  from  the  electrically  scanned  microwave 
radiometer  (ESMR)  aboard  Nimbus-6.  He  found  excellent  agreement  between  the 
observations  and  his  model’s  predictions  concerning  the  roughness  effect  but 
considerably  more  spread  in  the  foam-effect  results. 

Goodberlet  et  al.  (1989)  used  the  “D-matrix”  approach  within  their  algorithm  for 
retrieving  WS  measurements.  This  approach  attempts  to  find  a  linear  relationship 
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between  observed  brightness  temperatures  and  geophysical  parameters.  As  is  the  case 
for  retrieving  any  of  the  environmental  parameters,  any  rain  in  the  FOV  of  the 
radiometer  must  be  taken  into  account.  They  developed  a  rain  flag  criterion  for  their 
algorithm  that  is  sensitive  to  rain  rates  below  1  mm  hr'*  (rain  flag  “1”).  As  the  rain 
intensity  increases,  rain  flags  2  and  3  appear  to  reveal  areas  of  heavy  precipitation.  It 
follows  that  as  the  rain  rate  increases,  the  accuracy  of  their  retrievals  decreases. 

The  determination  of  cloud  liquid  water  over  the  oceans  using  microwave 
measurements  has  been  largely  neglected  (Bauer  and  Schluessel  1993;  Wilheit, 
personal  communication).  Alishouse  et  al.  (1990)  developed  a  purely  statistical 
method  to  determine  CLW  and  found  no  significant  correlation  between  observed 
SSM/I  brightness  temperatures  and  the  surface-based  CLW  calculations  when 
retrievals  were  performed  over  land  or  snow.  However,  they  found  a  high  correlation 
for  ocean-based  retrievals.  Greenwald  et  al.  (1993)  employed  a  simple  physical 
model  which  used  the  19  and  37-GHz  channels  of  the  SSM/I  to  retrieve  CLW 
estimates  over  the  oceans.  They  compared  their  physically  based  model  to  the  model 
developed  by  Alishouse  et  al.  They  found  that  the  Alishouse  et  al.  model 
underestimated  the  liquid  water  path  by  roughly  a  factor  of  2.  They  agreed  that  the 
statistical  model  detected  a  CLW  signal  but  the  underestimation  was  probably  due  to 
the  inability  of  that  model  to  completely  separate  the  influence  of  water  vapor  from 
the  effects  of  clouds. 

Optically  thick  paths  in  the  atmosphere  (containing  dense  clouds  or  rain) ,  i.e.,  the 
emission  and  scattering  regimes,  within  the  field  of  view  (FOV)  of  the  radiometer  can 
significantly  reduce  the  polarization  difference.  Additionally,  the  emission  of 
microwaves  from  cloud  droplets  and  rain  will  translate  as  an  increase  in  brightness 
temperatures  in  both  polarizations. 

Passive  microwave  measurements  reveal  information  about  rain  and  precipitation¬ 
sized  ice  in  the  FOV  of  the  SSM/I  sensors.  Prabhakara  et  al.  (1992)  used  the  37-GHz 
channels  to  compute  rain-rate  statistics  over  the  oceans.  They  found  that  at  this 
frequency,  where  the  FOV  at  the  Earth’s  surface  is  approximately  30  km,  rain  is  a 
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fairly  strong  emitter  of  microwave  radiation  while  scattering  by  ice  has  a  much 
smaller  effect.  At  frequencies  below  37  GHz,  where  the  FOV  is  larger  and  the 
volume  extinction  coefficient  is  weaker,  any  additional  information  about  rain 
appears  somewhat  redundant.  However,  at  85  GHz,  where  the  FOV  is  approximately 
15  km  and  the  extinction  is  strong,  direct  information  pertaining  to  rain  is  masked  by 
the  clouds.  At  this  frequency,  scattering  by  ice  becomes  very  noticeable. 

Nelkin  (1992)  examined  the  distribution  of  observed  brightness  temperatures  as 
measured  by  the  SSM/I  at  85  GHz  during  August  1987.  His  study  focused  on  several 
land  and  oceanic  regions  each  encompassing  an  area  of  5°  latitude  by  5°  longitude. 
For  each  area  for  the  entire  month,  he  examined  the  pairs  of  brightness  temperatures 
(vertical  and  horizontal),  calculated  each  pair’s  average  temperature,  (Tbh  +  Tbv)/2, 
along  with  their  polarization  difference  (Tbv  -Tbh  or  ATb)  and  plotted  these 
temperatures  on  two-dimensional  histograms  or  scatterplots.  The  observations  were 
plotted  on  a  Cartesian  grid  with  the  number  of  occurrences  of  each  observation 
indicated  by  a  shading  scale.  The  plots  were  displayed  with  the  average  temperatures 
on  the  vertical  axis  and  the  polarization  differences  on  the  horizontal  axis.  Figure  1  is 
an  example  of  the  resulting  plots  which  shows  two  ‘branches’.  The  vertically  oriented 
branch  is  indicative  of  the  scattering  regime  while  the  negatively  sloped  branch  is 
representative  of  the  absorption  regime. 

Next,  Nelkin  defined  a  somewhat  arbitrary  point  to  delineate  the  scattering  and 
absorption  regimes.  He  identified  a  point  at  the  top  of  the  “notch”  between  the  two 
branches  where  no  observations  were  recorded.  All  observations  to  the  right  of 
this  point  fell  into  the  absorption  branch;  all  observations  to  the  left  and  below  this 
point  were  classified  as  the  scattering  branch.  This  still  left  quite  a  few  observations 
unclassified  due  to  the  ambiguities  associated  with  the  source  of  the  atmosphere’s 
increased  opacity  (i.e.  scattering  from  raindrops  and  ice  particles  or  absorption  from 
clouds  and  water  vapor). 
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Nelkin  used  those  pairs  of  brightness  temperatures  in  the  absorption  branch  to 
restrict  retrievals  of  various  geophysical  parameters  using  the  technique  of  Wilheit  and 
Chang  (1980).  New  “restrictive”  algorithms  resulted  which  improved  the  accuracy  of 
the  retrieved  parameters.  However,  only  those  pairs  of  Tb  85  GHz  observations  and 
their  associated  lower  channel  observations  which  were  clearly  identified  as  belonging 
to  the  absorption  branch  could  be  input  into  these  new  algorithms.  Due  to  the 
increased  opacity  (from  rain  and  dense  clouds)  associated  with  many  weather  regimes, 
the  applicability  for  using  these  new  algorithms  was  severely  limited  (as  little  as  8% 
of  the  time  for  parts  of  the  ITCZ). 

4.  Research  objectives 

The  first  objective  of  this  investigation  is  to  determine  how  often  Nelkin’s 
improved  algorithms  can  be  applied  on  a  global  scale  after  classifying  all  observations 
as  belonging  to  either  the  absorption  or  scattering  regime.  These  algorithms  will 
constitute  a  simulation  study  for  the  performance  of  the  MIMR  which  is  tentatively 
scheduled  for  launch  in  the  year  2000.  The  result  of  my  first  objective  is  a  series  of 
global  maps  depicting  various  sensitivity  distributions  relating  the  polarization 
differences  between  the  85  GHz  channels  (Tbssv  -Tbssh)  with  the  opacity  of  the 
atmosphere  for  a  given  region. 

A  second  objective  is  to  understand  and  explain  the  sensitivity  distributions  given 
the  information  contained  in  the  CLW  distributions  over  selected  oceanic  regions. 
Few  studies  have  been  conducted  explaining  the  distribution  of  CLW  in  areas 
containing  nonprecipitating  clouds. 
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CHAPTER  II 
DATA 

The  data  used  in  this  study  consist  of  earth  located  brightness  temperatures  as 
measured  by  the  SSM/I  and  European  Center  for  Medium-range  Weather  Forecasts 
(ECMWF)  sea  surface  temperatures  generated  from  the  National  Center  for 
Atmospheric  Research  (NCAR)  data. 

The  SSM/I  data  cover  seven  months  (August  -  November  1987  and  January  - 

March  1993).  The  data  sets  are  written  in  a  binary  format  on  CD-ROMs  and  are 
arranged  in  either  hourly  or  daily  files.  The  files  contain  brightness  temperatures  for 
all  seven  channels,  time,  and  the  associated  location  of  an  observation  (latitude  and 
longitude).  In  order  to  keep  the  amount  of  data  to  a  manageable  level,  observation 
locations  are  reported  only  with  those  associated  with  the  low  frequency  channels  (i.e. 
A-scans  only,  64  per  scan).  Tg  85  observations  which  are  not  collocated  with  low 
frequency  channels  are,  for  our  purposes,  assigned  locations  to  an  adjacent  low 
frequency  channel  location  (see  Chapter  V  for  further  explanation  and  a  schematic 
representation). 

The  geographic  area,  for  this  study,  is  restricted  to  brightness  temperature 
measurements  over  the  oceans  between  60®  N  and  S  latitudes.  A  land  mask,  with  a 
grid  resolution  of  0.5°,  is  applied  to  enforce  this  restriction. 

Brighmess  temperatures  are  extracted  off  the  CD-ROMs  and  placed  into  monthly 
files.  These  files  are  further  subdivided  into  5°  x  5°  latitude/longitude  boxes.  This 
process  yields  1728  boxes  of  data  for  each  month  (land  and  ocean  combined). 

Boxes  that  contain  all  land  are  discarded  while  water  boxes  are  kept  for  further 
analysis.  Boxes  that  comprise  a  mixture  of  land  and  water  (i.e.  coastal  and  island 
areas)  are  classified  oceanic  the  water  area  is  >  10%  of  the  total  area.  Based  on  the 
above  criteria,  the  resulting  number  of  oceanic  boxes  is  1270. 
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Even  though  a  box  may  contain  as  little  as  10%  water,  only  those  observations 
that  fall  outside  the  land  mask  are  used.  This  ensures  land  contamination  is 
minimized.  However,  oceanic  observations  located  within  approximately  50  km  of 
land  are  subject  to  land  contamination  since  the  criteria  for  determining  water  vs.  land 
pixels  used  in  constructing  the  land  mask  is  imperfect. 

The  ECMWF  data  covers  the  time  period  from  mid-January  through  March  1993. 
Similar  to  the  SSM/I  data  sets,  ECMWF  surface  temperatures  are  also  written  in  a 
binary  format  and  stored  on  CD-ROMs.  Individual  files  contain  one  day  of  data. 
Observation  points  have  a  separation  of  2.5°  latitude  and  longitude  and  are  reported 
twice  daily  (OOZ  and  12Z). 
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CHAPTER  III 

ABSORPTION  VS.  SCATTERING 
OVER  THE  OCEANS 


1.  The  absorption  branch 

Wilheit  et  al.  (1994)  stated  that  brightness  temperatures  that  are  low  compared  to 
the  surface  temperature  and  the  associated  layers  of  the  atmosphere  through  which  the 
radiation  passes  are  produced  in  two  ways.  The  first  is  scattering  at  short  wavelengths 
(frequencies  >  60  GHz)  by  ice  particles  with  densities  and  sizes  comparable  to  rain. 
This  applies  particularly  well  to  the  SSM/I’s  85  GHz  channels.  The  second  is  by 
reflection  from  a  surface  that  has  a  high  reflectivity.  Such  is  the  case  with  an  ocean 
background  as  viewed  by  the  SSM/I  (Earth  incidence  angle  of  53.1°).  Since  the  view 
angle  is  far  from  nadir,  the  Fresnel  relations  from  electromagnetic  theory  (Griffiths 
1989)  require  that  the  reflectivity  for  horizontally  polarized  radiation  always  exceed 
the  reflectivity  for  vertically  polarized  radiation  (Fig.  2).  Conversely,  the  emissivity 
for  horizontally  polarized  radiation  will  always  be  less  than  the  vertical.  Therefore, 
the  vertically  polarized  brightness  temperature  will  always  be  greater  than  the 
horizontally  polarized  brightness  temperature. 

The  atmosphere  is  fairly  transparent  to  microwave  radiation  when  no  clouds  and 
little  or  no  water  vapor  is  present  (mild  absorption  case).  In  this  instance,  the 
brightness  temperatures  are  low  due  to  the  highly  reflective  property  of  the  ocean 
surface.  Referring  to  the  Fresnel  relations,  the  polarization  differences  (ATb  )  will  be 
large. 

As  the  opacity  of  the  atmosphere  increases,  through  the  addition  of  water  vapor 
and  clouds,  a  comparative  increase  in  the  average  brightness  temperatures  takes  place. 
The  absorption  and  concomitant  emission  of  microwave  radiation,  from  the  now 
moderate  to  strongly  absorbing  layers  of  the  atmosphere,  result  in  increased  brightness 


Incidence  Angle  (Degrees) 


FIG.  2.  Ocean  surface  reflectivity  as  a  function  of  incidence  angle  and 
polarization  for  85  GHz  and  surface  temperature  of  288  K.  (Ocean  surface  computed 
for  a  complex  index  of  refraction  of  3.23  +  1.88i). 
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temperatures.  Also,  the  ATb ’s  decrease  as  the  polarization  dependence  decreases  due 
to  the  increasing  opacity  of  the  intervening  atmosphere. 

Increasing  atmospheric  opacity  moves  the  Tb’s  along  a  nearly  straight  line  in  both 
polarizations  (from  lower  to  higher  average  Tb’s  and  concurrently,  larger  to  smaller 
ATb’s).  The  observations  fall  within  a  Tbv  ,  Tbh  space  (or  any  linear  transformation 
thereof)  which  can  be  depicted  as  the  negatively-sloped  absorption  branch  using  the 
average-difference  plots  (Fig.  1).  The  characteristic  slope  of  the  absorption  branch  is 
always  negative.  A  further  physical  argument  is  presented  in  section  4. 

2.  The  scattering  branch 

The  scattering  regime  can  be  characterized  by  small  differences  in  the  ATb’s  ,  and 
simultaneously,  lower  average  brightness  temperatures.  The  scattering  regime 
appears  as  a  nearly  vertically  sloped  branch  when  displayed  on  the  average-difference 
histograms  (Fig.  1). 

Low  brightness  temperatures  which  are  generated  by  scattering  will  depend  only 
slightly,  if  at  all,  on  the  polarization.  This  is  the  reason  for  very  small  ATb ’s  in  the 
scattering  branch.  It  is  the  observations  within  this  branch  that  must  be  eliminated  in 
order  to  obtain  reliable  algorithm  outputs.  Any  rain  in  the  retrieval  will  invalidate  the 
results. 

Following  the  technique  of  Nelkin  (1992),  average-difference  2-D  histograms 
are  generated  for  each  of  the  1270  oceanic  boxes  for  all  seven  months.  The  plots  (Fig. 
1  is  an  example)  graphically  depict  the  absorption  and  scattering  branches.  However, 
Nelkin’ s  notch  criterion  for  determining  the  borderline  between  both  branches  does 
not  account  for  a  large  number  of  observations  which  can  be  classified  as  belonging  to 
either  regime.  In  order  to  maximize  the  percentage  of  time  that  4-parameter  (CLW, 
PW,  SST,  and  WS)  algorithms  can  be  applied,  it  is  necessary  to  statistically  classify  as 
many  pair  of  85  GHz  (vertical  and  horizontal)  observations  as  possible  as  belonging 
to  either  the  absorption  or  scattering  regime.  This  entails  finding  a  method  to 
separate  the  observations  into  their  respective  regimes. 
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The  clue  to  solving  this  problem  lies  with  the  characteristic  slopes  associated  with 
each  branch.  The  slope  of  the  absorption  branch  is  always  negative  whereas  the  slope 
of  the  scattering  branch  is  vertical. 

The  problem  is  where  the  two  branches  join.  Here,  the  cloud  drops,  small 
raindrops,  and  water  vapor  are  all  contributing  to  the  increasing  brightness 
temperatures  via  their  absorption  and  emission  of  the  microwave  radiation.  The 
raindrops  have  yet  to  reach  a  size  sufficient  to  produce  noticeable  scattering. 

3.  Separation  of  the  branches 

The  following  objective  method  is  employed  to  delineate  the  observations  within 
each  branch.  First,  a  linear  regression  is  performed  using  the  method  of  least  squares 
to  find  the  slope  of  the  absorption  branch  for  a  given  5°  x  5°  oceanic  box. 
Specifically,  only  the  observations  which  can  be  unambiguously  identified  as 
belonging  to  the  absorption  branch  (i.e.  as  determined  by  Nelkin’s  notch  criterion)  are 
incorporated  into  the  regression.  The  standard  deviation  (a)  of  the  average  brightness 
temperatures  is  computed  concurrently.  Next,  the  resulting  regression  line  (which  runs 
through  the  center  of  the  absorption  branch)  is  extrapolated  to  cover  the  entire 
observed  dynamic  ATg  range  (0-72K). 

Before  proceeding  further,  an  assumption  is  made  that  the  average  brightness 
temperatures  within  the  absorption  branch  are  normally  distributed.  Boundary  lines 
are  then  displaced  ±  2a  from  the  regression  line.  Based  on  the  above  assumption, 
approximately  95%  of  all  observations  will  lie  between  the  boundary  lines.  All 
observations  that  fall  between  these  lines  are  classified  as  belonging  to  the  absorption 
regime.  Observations  outside  the  lines  are  classified  as  either  belonging  to  the 
scattering  regime  or  bad  data.  In  either  event,  for  the  purposes  of  this  study,  only  the 
observations  within  the  absorption  regime  are  relevant. 

The  histogram  in  Fig.  3  shows  the  result  of  applying  the  branch  separation 
method  to  the  data  used  in  plotting  Fig.  1.  This  box  lies  within  the  Intertropical 
Convergence  Zone  (TTCZ).  Approximately  17%  of  he  observations  fall  outside  the 
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absorption  branch.  Statistically,  this  is  equivalent  to  stating  that  it  rained 
approximately  17%  of  the  time  in  this  region  during  August  1987.  Climatologically, 
this  particular  statistic  appears  to  slightly  overestimate  (underestimate)  the  percentage 
of  scattering  (absorption)  related  brightness  temperature  observations.  This 
discrepancy  occurs  in  regions  which  commonly  experience  dense  clouds  and 
persistent  light  to  heavy  rainfall.  The  ITCZ,  South  Pacific  Convergence  Zone 
(SPCZ),  and  monsoon  regions  (i.e.  Bay  of  Bengal)  can  be  described  in  this  manner. 

However,  it  is  preferable  to  underestimate  (versus  overestimate)  the  number  of 
absorption  related  observations  since  the  main  drawback  is  a  reduction  in  the 
percentage  of  time  4-parameter  algorithms  can  be  applied.  An  overestimation  would 
allow  scattering  related  observations  into  the  algorithms.  The  errors  associated  with 
algorithm  outputs  would  increase.  However,  a  quantitative  assessment  of  these  errors 
has  yet  to  be  accomplished. 

Further  examples  (Figs.  4  and  5)  demonstrate  that  the  separation  method  works 
well  over  most  of  the  Earth’s  oceanic  regions.  Figures  4  and  5  are  representative  of 
the  midlatitudes  and  subtropical  high  belt  regions  respectively.  Two  distinct 
differences  between  the  midlatitudes  and  subtropics  versus  the  tropics  are  worth 
noting.  First,  the  plots  associated  with  the  midlatitudes  and  subtropics  reveal  a  higher 
percentage  of  observations  belonging  to  the  absorption  branch  than  do  those 
associated  with  the  tropics.  The  absorption  branch  in  the  midlatitude  box  (Fig.  4), 
accounts  for  98.5%  of  the  total  observations,  while  the  subtropical  box  (Fig.  5)  can  be 
described  as  a  single  branch  plot  (100%  absorption). 

Transient  extratropical  cyclones  and  their  associated  fronts  create  the  vertical 
updrafts  necessary  to  produce  convective  rain  and  ice  particles  which  are  associated 
with  the  small  percentage  of  scattering  observations  in  the  midlatitudes.  The  single 
branch  structure  of  the  subtropical  high  region  is  indicative  of  a  lack  of  vertical 
updrafts  necessary  to  produce  convective  rain.  Indeed, .  subtropical  highs  are  often 
characterized  as  large  areas  of  subsidence  which  inhibit  cloud  formation. 
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August  1987,  Midlatitude  (40-45'^N,  20-25®WI 
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FIG.  4.  SSM/I-observed  average  brightness  temperature  vs.  polarization  at  85 
GHz  for  August  1987  -  midlatitude  box.  (Data  between  diagonal  lines  constitutes  the 
absorption  branch  for  this  box). 
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August  1987,  Subtropical  High  (25-30‘'N,  20-25*’ VV) 
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FIG.  5.  SSM/I-observed  average  brightness  temperature  vs.  polarization  at  85 
GHz  for  August  1987  -  subtropical  high  box.  (Data  between  diagonal  lines  constitutes 
the  absorption  branch  for  this  box). 
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Second,  CLW  and  PW  vary  greatly  between  regions.  This  can  be  illustrated  by 
the  various  absorption  branch  lengths  and  widths.  Following  the  physical  description 
of  the  absorption  branch  in  section  2,  these  large  differences  of  CLW  and  PW  can  be 
explained  by  the  differing  climatic  regimes  associated  with  each  region.  The  ITCZ 
shows  the  least  variability  (short  and  narrow  branch),  whereas  the  midlatitude  box 
shows  the  greatest  variability  (long  and  wide  branch).  The  atmosphere  is  very 
opaque  along  the  ITCZ  due  to  the  extended  coverage  of  deep  convective  clouds.  The 
warming  of  the  brightness  temperatures  are  a  function  of  the  increasing 
opacity  of  the  intervening  atmospheric  layers.  Conversely,  variable  weather 
patterns  associated  with  the  midlatitude  box  (dynamic  frontal  boundaries  to  clear 
skies)  produce  a  wide  range  of  atmospheric  opacities. 

Even  though  the  branch  separation  method  works  well  over  the  majority  of 
the  Earth’s  oceans,  it  should  be  noted  that  sea  ice  (Fig.  6)  in  the  SSM/I’s  FOV 
yields  unusable  results.  It  is  apparent,  from  this  plot,  that  the  majority  of  observations 
form  a  vertical  sloping  branch.  It  should  also  be  evident  that  it  did  not  rain  99%  of 
the  time  during  February  1993  over  Hudson  Bay  (as  evidenced  from  the  scattering 
branch). 

Sea  ice  is  electrically  closer  to  land  (small  complex  index  of  refraction)  than 
water  which  gives  the  surface  a  low  reflectivity.  The  lowering  of  brightness 
temperatures  from  apparent  scattering  results  in  little,  if  any,  polarization  dependence 
and  thus  small  ATg  ’s.  Since  the  absorption  branch  carmot  be  detected,  any  areas 
containing  sea  ice  must  be  avoided  when  invoking  the  4-parameter  algorithms. 

This  also  suggests  that  small  ATb  ’s  are  a  sensitive  detector  for  the  presence  of  sea 
ice  and,  therefore,  can  be  used  to  filter  out  sea  ice  in  addition  to  rain  (reference 
Chapter  L  section  1). 

4,  Physical  justification  for  branch  separation  method 

Up  to  this  point,  a  statistical  argument  for  separating  the  absorption  and  scattering 
regimes  by  employing  Tg  85  observations  was  presented.  The  statistical  method  was 
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February  1993,  Hudson  Bay  (55-60°  N,  85-90°  W) 
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FIG.  6.  SSM/I-observed  average  brightness  temperature  vs.  polarization  at  85 
GHz  for  February  1993  -  Hudson  Bay  region.  (Branch  separation  method  not  effective 
due  to  presence  of  sea  ice). 
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FIG.  7.  Two-component  (single  layer  atmosphere  and  surface)  radiative  transfer 
model  (a)  and  effects  on  brightness  temperatures  as  a  result  of  varying  optical  depth 
(x)  (b).  Surface  emissivity  values  taken  from  figure  2.  Optical  depth  varied  from 
totally  opaque  (x=qo)  to  clear  sky  (x=0). 
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used  because  of  its  relative  simplicity  and  expeditious  processing  of  vast  quantities  of 
data  required  to  obtain  global  statistics.  However,  this  method  lacks  physical 
meaning.  Therefore,  a  simple  radiative  transfer  model  (RTM)  based  on  the  ERT  is 
used  to  show  that  the  statistical  method  is  physically  reasonable. 

A  simple  two-component  (single  layer  atmosphere  and  surface)  RTM  (Fig.  7a) 
can  be  constructed  to  illustrate  radiative  transfer.  This  model  does  not  account  for 
scattering  and  assumes  atmospheric  homogeneity.  For  a  given  surface  temperature 
(Tj),  view  angle  (0),  and  wavelength  (^),  by  varying  the  optical  depth  (x)  from  zero 
(no  atmosphere)  to  infinity  (completely  opaque),  brightness  temperatures  will  vary 
nearly  linearly  (in  both  polarizations)  from  some  computed  value  (no  atmosphere 
case)  to  the  cloud  thermometric  temperature  (totally  opaque  case)  (Fig.  7b). 

For  the  ITCZ  case,  this  temperature  range  is  indicated  by  the  lower  dotted  line  in 
Fig.  8  which  assumes  that  the  coldest  nonprecipitating  (non  scattering)  cloud  layer  to 
be  272  K.  If  the  cloud  temperature  is  assumed  to  be  close  to  the  surface 
temperature,  then  the  warmest  expected  temperature  (in  the  opaque  case)  would  be 
equivalent  to  that  of  the  surface.  This  temperature  range  is  indicated  by  the  upper 
dotted  line  in  Fig.  8. 

The  dotted  lines  in  Fig.  8  trace  out  a  triangle  with  the  zero  ATb  line  its  base.  The 
area  inside  this  triangle  represents  the  range  of  possible  values  which  could  occur 
by  allowing  attenuation  to  vary  from  zero  (no  atmosphere)  to  infinity  (totally  opaque) 
at  any  temperature.  The  underlying  physics  can  be  applied  globally,  however,  the 
expected  range  of  temperatures  (assuming  absorption  is  taking  place)  between  the 
surface  temperature  and  cloud  temperature  varies  latitudinally. 

The  temperature  range  is  dynamic  due  to  differing  cloud  types,  heights,  and 
thickness  which  are  functions  of  the  cloud’s  microphysics.  Cloud  microphysics  is 
associated  with  various  weather  regimes  which  in  turn  can  be  categorized  by  latitude. 

For  example,  in  the  tropics,  clouds  are  generally  convective  (large  vertical  extent) 
and  are  situated  globally  in  the  warmest  atmosphere.  This  condition  persists  to  the 
point  where  there  exists  relatively  little  variability  of  cloud  types  and  amounts. 
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FIG.  8.  Intercomparison  of  the  physical  model  and  the  statistical  method  to 
identify  the  absorption  branch  of  a  February  1993,  ITCZ  box.  (Data  within  dotted 
lines  associated  with  physical  model,  data  within  solid  diagonal  line  associated  with 
statistical  method). 


FIG.  9.  Intercomparison  of  the  physical  model  and  the  statistical  method  to 
identify  the  absorption  branch  of  a  February  1993,  midlatitude  box.  (Data  within 
dotted  lines  associated  with  physical  model,  data  within  solid  diagonal  line  associated 
with  statistical  method). 
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Therefore  a  dynamic  range  of  30  K  is  sufficient  to  account  for  this  variability  and 
serve  as  the  boundary  for  the  absorption  regime.  The  observations  lying  within  the 
model’s  triangular  boundaries  are  nearly  identical  to  those  which  are  defined  by  the 
statistical  method  (between  the  solid,  parallel  lines).  Since  the  areas  traced  out  by 
both  methods  are  similar,  this  serves  as  an  argument  that  the  statistical  method  is 
physically  reasonable. 

Another  example  which  illustrates  how  the  temperature  range  is  a  function  of 
cloud  physics  and  variability  is  the  midlatitudes.  The  midlatitudes  experience 
the  largest  variabilities  of  CLW  and  PW  due  to  transient  extratropical  cyclones. 
Large  daily  swings  in  atmospheric  opacity  (from  clear  skies  through  dense  overcasts) 
are  not  uncommon.  This  necessitates  a  larger  dynamic  temperature  range  of  45  K  to 
account  for  the  majority  of  possible  weather  regimes  (Fig.  9).  As  with  the  tropical 
region,  the  RTM  output  (area  within  the  dotted  triangular  region)  is  in  close 
agreement  with  the  statistical  method’s  output  (area  between  solid,  parallel  lines). 

As  long  as  a  sufficient  dynamic  temperature  range  is  specified,  the  RTM  output 
closely  mirrors  the  statistical  output.  A  reasonable  representation  for  comparing  the 
physical  and  statistical  methods  can  be  obtained  by  applying  the  dynamic  temperature 
range  criteria  found  in  Table  2.  This  table  differentiates  the  variability  of  atmospheric 
opacity  by  latitude. 
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TABLE  2.  Dynamic  temperature  range  necessary  for  application  of  physical 
model  to  statistics.  Climate  zones  based  on  analysis  of  February  1993  SSM/I  Tg  85 
GHz  observations  over  the  North  Pacific. 


Climate  Zones 

Latitude  Range 

(Degrees) 

Dynamic  Temperature 

Range  (K) 

High  Latitudes 

50-60 

40 

Midlatitudes 

30-50 

45 

Subtropics 

20-30 

40 

Tropics 

0-20 

30 

30 


CHAPTER  IV 

GLOBAL  SENSITIVITY  STUDY 

1.  Sensitivity  percentile  thresholds 

Atmospheric  opacity,  in  this  study,  is  described  in  terms  of  85  GHz  ATb’s.  As 
previously  mentioned,  increasing  atmospheric  opacity  yields  decreasing  ATb’s.  This 
allows  algorithm  sensitivities  to  be  described  in  terms  of  the  ATb  range. 

Regions  within  the  radiometer’s  FOV  where  the  atmosphere  is  relatively 
transparent  (i.e.,  ATb  ^  25  K)  can  be  identified  so  that  an  algorithm  can  be  made 
which  is  very  sensitive  to  a  specific  parameter  (CLW,  PW,  SST,  or  WS),  thereby 
eliminating  much  of  the  retrieval’s  uncertainty.  As  opacity  increases,  the  algorithm 
must  be  less  sensitive  to  the  desired  parameter.  This  results  in  a  larger  retrieval 
uncertainty. 

Four  percentile  thresholds  (10,  50,  75,  and  90)  have  been  chosen  to  represent  how 
often  algorithms  can  be  used  based  on  their  sensitivity.  The  method  to  determine  the 
distribution  of  ATb  sensitivities  involves  the  following  procedure  which  is  applied  to 
each  5°  x  5°  box  separately. 

First,  the  total  number  of  Tb  85  observations  are  tabulated.  Next,  only  the 
observations  within  the  absorption  branch  are  counted,  in  descending  order  of  ATb  s 
(from  right  to  left  along  the  horizontal  axis  in  Figs.  1  and  3-6).  When  the  number  of 
observations  within  the  absorption  branch  reaches  a  percentile  threshold  (percent  of 
the  total  observations),  the  corresponding  ATb  is  noted  and  assigned  as  the  ‘sensitivity 
value’  for  that  box  at  that  percentile  threshold. 

For  example,  referring  to  the  plot  in  Fig.  3,  the  total  number  of  observations  sums 
to  40,000.  Counting  observations  within  the  absorption  branch  (from  right  to  left), 
4,000  observations  (10  %  of  the  total)  have  ATb  s  greater  than  or  equal  to  19  K. 
Continuing  the  process,  20,000  observations  (50  %  of  the  total)  have  ATb ’s  greater 
than  or  equal  to  14  K.  This  process  continues  until  all  four  percentile  thresholds  are 
calculated.  If  the  number  of  observations  within  the  absorption  branch  is  less  than 
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what  is  required  for  a  particular  threshold  (i.e.  there  is  fewer  than  36,000  observations 
to  attain  90  %  in  Fig.  3),  then  a  default  ATg  value  of -1  K  is  assigned. 

2.  Global  sensitivity  distributions 

After  sensitivity  values  are  assigned  to  all  1270  oceanic  boxes,  global  sensitivity 
maps  for  each  of  the  four  percentile  thresholds  are  constructed  (Figs.  10-13).  Each 
figure  shows  a  monthly  (August  1987)  composite,  for  a  given  percentile  threshold,  of 
the  1270  regions  and  their  assigned  sensitivity  values. 

Based  on  the  assigned  shading  scale,  sensitivity  values  (or  equivalently,  85  GHz 
polarization  differences)  indicate  relatively  optically  transparent  (blue,  >25  K)  to 
relatively  opaque  (red,  0-5  K)  regions.  The  percentiles  indicate  how  often  algorithms 
can  be  invoked.  This  depends  on  their  individual  sensitivity  for  retrieving  a  specific 
geophysical  parameter.  For  example,  in  Fig.  10,  a  relatively  sensitive  algorithm 
(sensitivity  value  >20  K)  can  be  employed  at  least  10%  of  the  time  over  the  majority 
of  the  Earth’s  oceans.  The  only  exceptions  appear  to  be  along  portions  of  the  ITCZ 
and  monsoon  regions  (Bay  of  Bengal  and  Arabian  Sea)  where  dense  clouds,  water 
vapor,  ice  particles,  and  persistent  rain  render  the  atmosphere  sufficiently  opaque  to 
lower  sensitivity  values  of  15-20  K.  The  same  logic  applies  for  the  remaining  three 
percentile  thresholds. 

A  by-product  of  the  global  maps  is  a  strong  correlation  between  the  ATg’s  and 
large  scale  weather  features.  This  is  a  good  cross  check  for  verifying  the  reliability  of 
the  entire  global  sensitivity  mapping  process  since  clouds,  ice  particles,  water  vapor, 
and  rain  are  the  main  ingredients  for  triggering  a  change  in  atmospheric  opacity  at  85 
GHz.  At  the  50th  percentile  (Fig.  11),  moderately  opaque  (sensitivity  values  15-25 
K)  regions  appear  in  the  extreme  South  Atlantic  Ocean,  while  the  majority  of 
all  southern  hemispheric  oceans  remain  relatively  transparent  (>  25  K). 

At  first  glance,  this  might  suggest  that  the  increased  opacity  is  a  result  of 
increased  cloudiness,  but  inspection  of  the  individual  plots,  of  these  relatively 
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FIG.  10.  Global  sensitivity  map  depicting  85  GHz  polarization  difference  values 
for  August  1987  at  the  10th  percentile  threshold.  (Shading  scale  applies  only  to  the 
1270  oceanic  regions;  land  regions  default  to  white). 
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FIG.  11.  Global  sensitivity  map  depicting  85  GHz  polarization  difference  values 
for  August  1987  at  the  50th  percentile  threshold.  (Shading  scale  applies  only  to  the 
1270  oceanic  regions;  land  regions  default  to  white). 
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FIG.  12.  Global  sensitivity  map  depicting  85  GHz  polarization  difference  values 
for  August  1987  at  the  75th  percentile  threshold.  (Shading  scale  applies  only  to  the 
1270  oceanic  regions;  land  regions  default  to  white). 
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FIG.  13.  Global  sensitivity  map  depicting  85  GHz  polarization  difference  values 
for  August  1987  at  the  90th  percentile  threshold.  (Shading  scale  applies  only  to  the 
1270  oceanic  regions;  land  regions  default  to  white). 
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opaque  regions,  reveal  similar  profiles  to  those  seen  in  Fig.  6.  Therefore,  based 
on  the  85-90  GHz  polarization  differences,  these  regions  are  further  examples  of 
SSM/I  observations  contaminated  by  sea  ice,  and  as  such,  cannot  be  incorporated 
into  the  global  sensitivity  distributions. 

At  the  75th  percentile  (Fig.  12),  portions  of  the  Bay  of  Bengal  are  extremely 
opaque  (sensitivity  values  <  0-5  K).  Specifically,  a  numerical  analysis  of  the  box 
which  indicates  a  value  of -1  K  (10  -  15°N,  95-100°E)  reveals  73.8%  absorption  based 
observations.  This  statistic  is  equivalent  to  stating  that  it  rained  26.2%  of  the  time  in 
that  region  during  August  1987.  To  keep  this  in  perspective,  it  would  have  to  rain  an 
average  of  six  hours  per  day  to  achieve  this  absorption/scattering  ratio.  Since  it  is 
during  the  monsoon,  this  statistic  seems  plausible  if  not  completely  accurate. 

Several  additional  meteorological  features  appear  in  fig.  12.  First,  the  SPCZ  is 
visible,  via  small  polarization  differences,  from  the  eastern  equatorial  Pacific  through 
the  South  Central  Pacific.  Next,  the  majority  of  the  eastern  North  Pacific  and  Atlantic 
Oceans  also  reveal  rapidly  decreasing  ATg’s.  In  these  regions,  the  increasing 
atmospheric  opacity  is  most  likely  associated  with  transient  midlatitude  cyclones 
which  feed  off  the  warm  waters  carried  by  the  Kuroshio  and  Gulf  Stream  currents. 

Finally,  there  are  many  boxes  which  remain  relatively  transparent  (  >20  K)  at  the 
75th  percentile  threshold  to  include  the  southern  hemispheric  oceans  between  20  and 
60°  S.  Transparent  areas  exist  along  the  western  coasts  of  North  America  and  Africa 
(minus  the  ITCZ).  It  is  in  these  areas  where  subsidence,  associated  with  subtropical 
highs,  is  a  major  determining  factor  for  the  scarcity  of  clouds  and  precipitation. 

The  global  sensitivity  map  of  the  90th  percentile  (Fig.  13)  continues,  as  expected, 
to  show  a  decrease  of  polarization  differences  through  all  boxes.  It  is  only  in  the  most 
opaque  regions  (black  (-1  K)  and  red  (0-5  K))  where  no  algorithms  can  be  expected  to 
produce  reliable  results  a  majority  of  the  time.  Specifically,  the  areas  of  concern 
include  the  TTCZ,  monsoon  regions,  portions  of  the  SPCZ,  and  extreme  South 
Atlantic  Ocean  (due  to  sea  ice  formation).  However,  this  map  shows  that  the  least 
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sensitive  algorithms  (those  which  function  down  to  a  A  Tg  of  5  K)  can  be  utilized  up 
to  90%  of  the  time  over  the  majority  of  the  oceans. 

An  alternative  method  for  displaying  the  global  sensitivity  distributions  involves 
focusing  on  the  latitudinal  variability  of  polarization  differences.  The  graphs  in  Figs. 
14  and  15  contain  a  plethora  of  information  relating  to  both  inter  and  intralatitudinal 
variability.  Each  graph  displays  statistics  (maximum  and  minimum  ATb’s,  mean  ATg, 
and  +  one  standard  deviation)  for  24  separate  sets  of  data. 

Each  data  set  contains  oceanic  boxes  grouped  along  latitudinal  bands  in 
increments  of  5°  from  60°  N  to  S.  Since  boxes  have  a  width  of  5°  longitude,  a 
maximum  of  72  boxes  is  possible  within  a  5°  latitudinal  band  around  the  Earth.  The 
actual  number  of  boxes  used  in  computing  the  statistics  for  a  given  latitudinal  band 
varies  and  can  be  considerably  less  especially  in  the  Northern  Hemisphere  where  the 
majority  of  Earth’s  land  masses  reside. 

At  the  10th  and  50th  percentile  thresholds,  there  exist  a  symmetric  pattern  to  the 
data  between  hemispheres.  This  pattern  implies  an  increase  in  atmospheric  opacity 
with  decreasing  latitude,  as  evidenced  from  the  mean  values.  The  summer  (Northern) 
Hemisphere  atmosphere  is  more  opaque  as  seen  from  its  lower  mean  ATg  values.  The 
intralatitudinal  variability  (±  a)  is  similar  for  both  hemispheres  with  one  notable 
exception;  Izirger  variability  between  30  and  40°  N  as  compared  to  other  latitudes. 
This  increased  variability  could  be  partially  due  to  the  extreme  differences  in  weather 
regimes  associated  with  the  subtropical  highs  (eastern  portion  of  ocean  basins)  and 
warm  ocean  currents  (western  portion  of  ocean  basins). 

At  the  75th  percentile  (Fig.  15),  the  symmetry  begins  to  break  down.  A  more 
complex  pattern  emerges  which  adds  information  to  the  rather  simple  hemispheric 
pattern  suggested  in  the  above  paragraph.  The  summer  hemisphere  continues  to  show 
a  higher  atmospheric  opacity;  however,  with  the  exception  of  the  tropics,  opacity  no 
longer  increases  as  latitude  decreases.  The  presence  of  the  subtropical  high  belt  can 
again  be  seen  by  its  large  ATg  variability.  The  global  ATg  variability  increases  for  all 
latitudes  from  the  75th  to  90th  percentiles  while  a  peak  mean  value  between  30  and 
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FIG.  14.  Statistics  displaying  the  latitudinal  variability  of  the  sensitivity 
distributions  for  August  1987  at  the  10th  (top)  and  50th  (bottom)  percentile 


thresholds. 
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FIG.  15.  Statistics  displaying  the  latitudinal  variability  of  the  sensitivity 
distributions  for  August  1987  at  the  75th  (top)  and  90th  (bottom)  percentile 


thresholds. 
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35°  N  (Fig.  15)  indicates  the  most  transparent  atmospheric  region  in  the  Northern 
Hemisphere. 

The  graph  displaying  the  90th  percentile  statistics  can  be  described  as  having  a 
skewed  distribution  (the  -a  exceeds  the  minimum  ATb  value)  for  the  boxes  in  the  data 
set  ranging  from  5  - 10°  N.  It  is  no  coincidence  that  this  happens  to  be  the  heart  of  the 
ITCZ. 

Two  factors  contributing  to  the  skewness  can  be  explained.  First,  the  minimum 
ATb  value  (•!  K)  is  an  artificial  lower  boundary  (refer  to  section  1,  this  chapter).  The 
majority  of  the  observations  within  the  ITCZ  are  extremely  opaque  and  warrant  the  -1 
K  value.  This  leads  to  an  extremely  small  mean  ATb  (1.7  K).  Second,  a  few  boxes 
are  associated  with  ATb’s  as  large  as  16  K.  This  computes  to  a  fairly  large  variability 
and  a  standard  deviation  of  ±  5.  Hence,  the  <s  value  (-3.3  K)  exceeds  the  artificial 
minimum  ATb  value  (-1  K). 

As  previously  mentioned,  a  good  cross  check  for  the  global  sensitivity  maps  are 
large  scale  weather  features.  In  addition  to  finding  latitudinal  variations  within  the 
sensitivity  distributions,  comparing  summer  and  winter  scenarios  yields  further 
information  regarding  seasonal  migration  of  large  scale  weather  patterns  (Figs  16-23). 
These  graphs  answer  the  questions;  “Do  seasonal  variations  in  weather  patterns  alter 
the  sensitivity  distributions?”  and  if  so,  “How  does  this  affect  the  applicability  of  our 
algorithms  for  retrieving  various  geophysical  parameters”  ? 

For  the  seasonal  comparison,  data  from  August  1987  and  February  1993  are 
analyzed  to  investigate  changes  that  occur  in  the  sensitivity  distributions.  At  the  10th 
percentile  (Fig.  16),  several  major  differences  are  evident.  First,  the  southward  shift 
of  the  ITCZ  (page  23,  Crowley  and  North  1991)  during  northern  winter  (January  - 
March)  is  readily  apparent.  Associated  with  the  ITCZ  shift,  the  monsoon  region  is 
experiencing  its  dry  season  (page  35,  Wallace  and  Hobbs  1977).  In  sharp  contrast  to 
the  summer  scenario,  the  skies  over  the  Arabian  Sea  and  Bay  of  Bengal  are  mostly 
cloud  free.  This  equates  to  large  polarization  differences  for  those  regions  which 
allows  the  use  of  algorithms  that  are  sensitive  to  CLW,  PW,  SST,  and  WS.  Also, 
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FIG.  20.  Statistics  comparing  the  latitudinal  variability  of  the  sensitivity 
distributions  for  August  1987  (top)  and  February  1993  (bottom)  at  the  10th  percentile 
threshold. 
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August  1987  (50th  percentile) 
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FIG.  21.  Statistics  comparing  the  latitudinal  variability  of  the  sensitivity 
distributions  for  August  1987  (top)  and  February  1993  (bottom)  at  the  50th  percentile 


threshold. 
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FIG.  22.  Statistics  comparing  the  latitudinal  variability  of  the  sensitivity 
distributions  for  August  1987  (top)  and  February  1993  (bottom)  at  the  75th  percentile 


threshold. 
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FIG.  23.  Statistics  comparing  the  latitudinal  variability  of  the  sensitivity 
distributions  for  August  1987  (top)  and  February  1993  (bottom)  at  the  90th  percentile 


threshold. 
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several  regions  in  the  Northern  Hemisphere  are  p>otentiaI  regions  of  sea  ice  formation: 
the  Sea  of  Okhotsk  (55  -  60°  N,  150  -  155°  E),  Hudson  Bay  (55  -  60°  N,  80  -  90° 
W),  and  the  Gulf  of  St.  Lawrence  (45  -  50°  N,  55  -  60°  W).  Examination  of  the 
individual  scatterplots,  such  as  Hudson  Bay’s  (Fig.  6),  reveal  the  characteristic 
scattering  branch  with  little  or  no  absorption  based  observations.  As  a  result,  none  of 
the  4-parameter  algorithms  are  expected  to  give  reliable  results  for  these  regions. 

At  the  50th  percentile  (Fig.  17),  the  polarization  differences  for  the  extreme  South 
Pacific  Ocean  in  February  exceed  25  K.  This  is  in  direct  contrast  to  the  August  plot 
which  depicts  regions  of  decreasing  ATg’s  in  the  South  Pacific  Ocean  between  55  and 
60°  S.  Sea  ice  does  not  extend  that  far  north  during  southern  summer  (January- 
March).  As  a  result,  the  more  sensitive  algorithms  can  be  used  in  these  regions  at 
least  50%  of  the  time. 

At  the  75th  percentile  (Fig.  18),  two  distinct  differences  emerge  between  the 
summer  and  winter  plots.  First,  the  polarization  differences  along  the  SPCZ  have 
decreased  significantly  in  February.  Interestingly,  this  decrease  exceeds  the  one 
shown  for  August.  The  rapidly  decreasing  ATb’s  are  associated  with  the  rapidly 
intensifying  SPCZ  during  southern  summer  (Vincent  1994).  Increased  convergence 
producing  denser  (convective)  clouds,  rain,  and  ice  particles  adds  to  the  atmospheric 
opacity  so  that  only  the  least  sensitive  algorithms  can  be  invoked. 

The  second  difference  involves  the  sharp  contrast  between  the  Northern  and 
Southern  Hemispheres.  The  regions  with  the  largest  polarization  differences  are 
associated  with  their  respective  winter  hemisphere.  Subsidence,  fewer  clouds,  and 
decreased  rainfall  are  associated  with  large  areas  of  high  pressure.  This  leads  to 
relatively  optically  transparent  atmospheric  conditions.  As  a  result,  the  more 
sensitive  algorithms  (>25  K)  are  usable  at  least  75%  of  the  time  over  the  winter 
hemisphere  oceans. 

As  with  all  the  global  maps,  the  trend  of  decreasing  ATb’s  continues  through  the 
90th  percentile  (Fig.  19).  The  most  striking  feature  for  comparison  involves  the 
SPCZ.  During  February,  the  convective  intensity  throughout  the  SPCZ  is  much 
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greater  than  August.  As  a  result,  even  the  least  sensitive  algorithms  (0-5  K)  are  not 
usable  in  the  SPCZ  during  southern  summer.  However,  at  least  90%  of  the  time,  the 
majority  of  the  remaining  global  oceanic  regions  are  usable  for  algorithms  with 
sensitivity  thresholds  above  5  K. 

Comparing  latitudinal  variability  between  seasons  (Figs.  20-23)  yields  some 
distinct  similarities  and  differences.  At  the  10th  percentile  (Fig.  20),  both  graphs 
show  strong  symmetry  between  hemispheres.  In  both  cases,  the  winter  hemisphere  is 
associated  with  larger  mean  ATb’s  .  Additionally,  both  seasons  show  a  rapid  decrease 
of  polarization  differences  with  decreasing  latitude.  The  graphs  depicting  the  seasonal 
variability  at  the  50th  percentile  (Fig.  21)  continue  to  show  strong  hemispheric 
symmetry  with  the  exception  of  the  region  north  of  40°N  during  February.  Here,  a 
reverse  of  the  decreasing  ATb’s  with  decreasing  latitude  is  noted. 

There  are  two  possible  explanations  for  this  discrepancy.  First,  several  regions 
of  sea  ice  are  present  (as  observed  by  the  presence  of  several  regions  of  extremely  low 
minimum  ATb  values).  The  arithmetic  mean  values  are  lower  than  expected  since  the 
mean  is  strongly  influence  by  these  outliers.  Another  potential  contributor  to  the 
reverse  ATb  trend  is  associated  with  the  weather  regime.  Sufficient  numbers  of 
transient  extratropical  cyclones  traversing  these  northern  latitudes  could  effectively 
lower  the  ATb’s  due  to  their  extensive  cloud  shields  and  associated  precipitation. 

The  overall  decrease  in  ATb’s  continues  throughout  the  75th  and  90th  percentile 
thresholds  (Figs.  20  and  21). 
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CHAPTER  V 

REGIONAL  SENSITIVITY  STUDY 

1.  Enlarging  the  Tg  85  GHz  FOV 

As  mentioned  in  Chapter  I,  the  SSM/I’s  and  MIMR’s  measurement  geometries  are 
such  that  there  are  many  high  resolution  Tg  85-90  observations  distributed  around  the 
low  resolution  channels.  Selecting  the  appropriate  algorithm  for  a  particular  parameter 
will  be  driven  by  the  lowest  value  ATg  observed  at  high  resolution  within  the  low 
channel’s  FOV. 

In  general,  the  smaller  the  ATg  value  the  more  robust  the  algorithm  needs  to  be. 
Unfortunately,  a  robust  algorithm  implies  a  relatively  large  uncertainty  in  the  retrievals 
compared  to  those  found  in  more  sensitive  algorithms.  Nelkin  (1992)  determined  the 
uncertainties  associated  with  the  4  -  parameter  algorithms  as  a  function  of  the 
minimum  85-90  GHz  ATg  (Table  3).  As  a  result  of  a  subsequent  error  found  with  his 
random  number  generator,  the  uncertainty  values  have  been  increased  by  a  factor  of 
1 .5  to  compensate  for  the  discrepancy. 

Since  we  are  taking  the  minimum  ATg  ,  if  there  is  any  variability  in  the 
observations,  the  sensitivity  distributions  will  shift  downward.  Also,  because  the 
extent  of  this  degradation  is  not  known,  a  procedure  to  simulate  decreasing  85  GHz 
channel  resolution  is  presented. 

This  procedure  involves  analyzing  progressively  larger  groups  of  collocated  Tg 
85  observations.  Groups  of  4  (2  along-scan  x  2  cross-scan),  16  (4  x  4),  and  64  (8  x  8) 
are  chosen  since  these  figures  are  multiples  of  the  number  of  total  observations  (128) 
sampled  along  a  scan  line.  Thus,  data  manipulations  and  calculations  are  performed 
without  the  need  to  incorporate  any  type  of  weighting  or  averaging  scheme  to  account 
for  any  ‘left  over’  observations.  It  is  of  interest  to  note  that  groups  of  4  and  16  have 
comparable  resolutions  with  the  37  and  22  GHz  channels  respectively.  The  group  of 
64  compares  favorably  to  the  resolution  that  a  10.65  GHz  channel  would  have. 
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TABLE  3.  4  -  parameter  algorithm  sensitivities  as  a  function  of  polarization 
difference.  Sensitivity  values  taken  from  Nelkin  (1992)  and  adjusted  to  compensate 
for  errors  in  his  random  number  generator. 


Min.  85-90 

GHz  ATb  (K) 

PW  (kg/m") 

CLW  (kg/m") 

SST  (K) 

WS  (m/s) 

0 

1.740 

.069 

1.357 

2.51 

18 

1.035 

.041 

1.152 

1.90 

20 

1.005 

,039 

1.153 

1.85 

22 

0.975 

.038 

1.122 

1.76 
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Once  observations  are  grouped,  the  pair  with  the  smallest  ATg  (from  a  group  of 
4,  16,  or  64)  is  identified.  This  pair  indicates  the  area  (85  GHz  FOV)  containing  the 
most  opaque  atmosphere  within  that  group.  The  remaining  pairs  in  the  group  are 
dropped  from  further  analysis,  while  the  smallest  ATg  value  is  retained  to  represent 
the  total  area  (FOV)  covered  by  the  group.  The  latitude  and  longitude  of  the  low 
channel  (A-scan)  observation  point  corresponding  most  closely  to  the  group’s 
geographic  center  is  then  assigned  this  ATg  value.  Based  on  the  assigned  latitude  and 
longitude,  the  ATg  value  is  then  placed  into  the  appropriate  5x5  degree  box.  From 
this  point  on,  the  data  are  processed  in  the  manner  described  in  Chapter  in. 

Figures  24,  25,  and  26  illustrate  the  scan  geometry  of  the  SSM/I.  However,  it 
should  be  noted  that  they  are  not  drawn  precisely  to  scale  and  that  the  actual  scans  are 
conical,  not  linear.  Large  X’s  denote  locations  where  the  radiometer  records 
observations  for  all  seven  channels.  Small  X’s  indicate  locations  where  only  the  85 
GHz  channels  are  sampled.  Small  ovals,  encircling  individual  observation  points, 
depict  the  approximate  spatial  resolution  of  the  85  GHz  channels  (15  km).  Large 
ovals  illustrate  the  approximate  37  and  22  GHz  FOV’s  (37  x  38  km  and  60  x  47  km. 
Figs.  24  and  25,  respectively).  Shaded  ovals  (dark  and  light)  designate  observation 
locations  representing  specific  groups  (4,  16,  or  64).  The  dark  shaded  ovals  indicate 
the  low  channel  observation  point  selected  to  represent  the  geographic  center  of  the 
group.  The  relative  horizontal  resolutions  of  the  19,  10,  and  6  GHz  channels  (69  x  50 
km,  120  X  85  km,  and  190  x  120  km,  respectively)  are  compared  to  each  other  and 
overlaid  onto  a  group  of  64  (Fig.  26). 

The  precise  center  of  a  group  does  not  directly  correspond  to  an  observation  point 
since  the  groups  are  composed  of  an  even  number  of  ‘rows’  and  ‘columns’  (cross-scan 
and  along-scan).  However,  aggregate  FOV’s  composed  of  various  85  GHz  groups 
compare  favorably  with  low  channel  FOV’s.  Figure  24b  shows  the  approximate  37 
GHz  FOV  and  its  adjoining  85  GHz  observations.  It  appears  that  as  many  as  nine  85 
GHz  FOV’s  (totally  or  partially)  are  incorporated  within  the  37’s  FOV.  However,  by 
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A  scan 
B  scan 
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FIG.  24.  Simplified  scan  geometry  of  the  SSM/I  with  relative  FOV’s  for  37  and 
85  GHz.  Group  of  4  Tb  85  FOV’s  highlighted  in  (a).  Tb 37  FOV  overlaid  onto 
several  Tb  85  FOV’s  to  show  relative  area  coverage  (b).  Tb  37  FOV  compared  to  a 
group  of  4  Tb  85  FOV’s  (c). 


55 


□1 

B 

^^9 

B 

Mo 

ma 

«!« 

SSI 

RHfili 

■ia^? 

S 

Aid 

«j« 

1 

1 

M 

w 

mm 

MW 

Ki 

iaKOHi 

Mfi! 

m 

FIG.  25.  Simplified  scan  geometry  of  the  SSM/I  with  relative  FOV’s  for  22  and 
85  GHz.  Group  of  16  Tb  85  FOV’s  compared  to  a  Tb  22  FOV  (left  panel). 
Tb22  FOV  overlaid  onto  several  Tb  85  FOV’s  showing  relative  area  coverage 
(right  panel). 


A  scan 


FIG.  26.  Simplified  scan  geometry  of  the  SSM/I  with  relative  FOV’s  for  6,  10 
19,  and  85  GHz.  Group  of  64  Tg  85  FOV’s  compared  to  Tb  19,  10,  and  6  FOV’s. 
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centering  the  37  GHz  FOV  on  a  group  of  four  (Fig.  24c),  it  is  apparent  that  the  total 
85  GHz  coverage  is  comparable  to  the  37’s  FOV.  Similarly,  as  many  as 
twenty-one  85  GHz  FOV’s  are  completely  or  partially  enclosed  within  the  22  GHz 
FOV  (Fig.  25b).  Once  again,  by  centering  the  22  GHz  FOV  on  a  group  of  16 
(Fig.  25a),  it  can  be  seen  that  a  group  of  16  has  an  equivalent  resolution  to  the  22 
GHz  channel. 

Finally,  the  19,  10,  and  6  GHz  FOV’s  are  overlaid  onto  a  group  of  64  (Fig.  26). 
The  area  covered  by  this  group  compares  most  closely  with  that  of  the  10  GHz  FOV, 
while  the  majority  of  the  6  GHz  channel  FOV  is  represented. 

Since  the  10  and  6  GHz  channels  are  not  supported  by  the  SSM/I,  their  respective 
resolutions  are  determined  from  the  existing  SSM/I  channels  and  then  scaled  based  on 
their  wavelengths.  Holding  all  other  factors  constant,  such  as  slant  range  (distance 
from  the  radiometer  to  the  ground  observation  point)  and  antenna  aperture  diameter, 
there  exists  a  nearly  linear  relationship  between  wavelength  and  resolution. 
Specifically,  resolution  increases  (smaller  FOV)  as  wavelength  decreases.  Although 
resolutions  will  increase  (by  a  factor  of  2  -  3  times)  aboard  the  MIMR,  the  channels 
relative  resolutions  will  remain  nearly  constant.  Therefore,  the  conclusions  resulting 
from  this  experiment  can  be  used  in  a  simulation  study  for  the  performance  of  the 
MIMR. 

2.  Effects  of  enlarging  the  FOV 

Since  the  37  GHz  resolution  (group  of  4)  is  coarser  than  the  85  GHz,  we  see  a 
degradation  of  sensitivities  (progressively  lower  ATg’s)  throughout  the  entire  globe 
(Figs.  27  and  28).  Similarly,  groups  of  16  yield  further  declines  in  the  global 
sensitivity  distributions  (Figs.  29  and  30),  while  groups  of  64  have  the  least  sensitive 
(most  opaque)  distributions  (Figs.  3 1  and  32). 

These  global  portrayals  of  the  85  GHz  sensitivity  distributions  are  primarily 
qualitative  in  the  sense  that  they  show  large-scale  atmospheric  weather  features  (in 
terms  of  atmospheric  opacity  or  optical  depth)  associated  with  both  cloud  mass  and 
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FIG.  27.  Global  sensitivity  maps  depicting  groups  of  4,  Tb  85  GHz  observations 
and  their  resulting  polarization  difference  values  at  the  10th  (top)  and  50th 
(bottom)  percentile  thresholds  for  August  1987. 


O  30E  60E 

90E 

120E 

ISOE  180  ISOM 

LONGITUDE 

120M 

■ 

I'iIvS-' 

Ilii 

ftugB7<V5th_»cit  i  le> 
4_observal ions 

<o 

o 

5  lO  15 

TbSSv  -  TbBShi 

20 

<K> 

>2! 

A:  li 


O  30E  60E  90E  ±20E  ISOE  180  150M  12014  SOM  60H  SON 

LONGITUDE 


ftug8V<90th_pct i le> 
4__of9£rerv-at  ions 


5  IQ  IS  20  >25 

TbSSv  -  TbBSh  <K) 


FIG.  28.  Global  sensitivity  maps  depicting  groups  of  4,  Tb  85  GHz  observations 
and  their  resulting  polarization  difference  values  at  the  75th  (top)  and  90th  (bottom) 
percentile  thresholds  for  August  1987. 
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FIG.  29.  Global  sensitivity  maps  depicting  groups  of  16,  Tg  85  GHz  observations 
and  their  resulting  polarization  difference  values  at  the  10th  (top)  and  50th  (bottom) 
percentile  thresholds  for  August  1987. 
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FIG.  32.  Global  sensitivity  maps  depicting  groups  of  64,  Tg  85  GHz  observations 
and  their  resulting  polarization  difference  values  at  the  75th  (top)  and  90th  (bottom) 
percentile  thresholds  for  August  1987. 
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water  vapor.  The  progressive  degradation  of  the  sensitivity  distributions  is  a  direct 
result  of  selecting  the  smallest  ATg  value  to  represent  increasingly  larger  groups  of  Tg 
85  observations  and  their  correspondingly  larger  FOV’s. 

Similar  to  the  cases  where  every  available  Tb  85  observation  is  utilized  to 
create  global  sensitivity  maps  (Figs.  16-19),  the  trend  of  decreasing  ATb’s  (from  the 
10th  through  90th  percentile  thresholds)  is  evident  no  matter  what  group  size  is 
chosen. 

This  is  due  to  the  maimer  in  which  ATb’s  are  computed  and  has  no  other 
significant  meaning,  physical  or  otherwise.  However,  the  trend  of  decreasing  ATb’s  , 
between  the  various  group  sizes  as  evidenced  from  the  global  profiles  in  figures  27- 
32,  has  physical  meaning.  To  some  extent,  these  figures  show  how  changing 
wavelengths  affects  resolution  which,  in  turn,  affects  the  sensitivity  distributions.  The 
exact  quantitative  measure  of  this  wavelength/resolution  effect  cannot  be  determined 
from  these  global  maps. 

Once  again,  an  alternative  method  for  displaying  the  results  can  be  shown  in 
terms  of  the  latitudinal  variability  of  the  polarization  differences  (Figs.  33-40).  This 
time,  the  attention  is  focused  on  comparing  the  intralatitudinal  variability  (±  ct  from 
the  mean)  between  the  various  Tb  85  group  sizes  at  specific  percentile  thresholds.  For 
example.  Figs.  33  and  34  depict  the  latitudinal  variability  of  the  polarization 
difference  at  the  10th  percentile.  Figure  33  contains  the  graphs  of  all  Tb  85 
observations  and  groups  of  four  during  August  1987.  Figure  34  contains  the  graphs 
for  groups  of  16  and  64  during  the  same  period.  Assuming  a  normal  or  near  normal 
distribution  of  ATb’s  along  a  given  latitude  belt,  an  interval  of  ±a  accounts  for 
approximately  68%  of  the  total  observations  (Ott  1993)  along  that  latitude.  This 
allows  for  a  more  quantitative  assessment  of  the  impact  that  resolution  plays  on  the 
sensitivity  distributions. 

Small  differences  exist  between  the  graphs  in  Fig.  33  that  represents  a  decrease  in 
resolution  by  approximately  a  factor  of  2  (when  considering  the  4-observation  graph). 
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August  1987  (10th  percentile,  all  observations) 
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August  1987  (10th  percentile,  4  observations) 
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FIG.  33.  Latitudinal  variability  of  the  Tb  85  GHz  polarization  differences  at  the 
lOth  percentile  threshold  for  August  1987  (all  Tb  85  observations  (top)  and  groups  of 
4  (bottom)).  Variability  is  ±  1  standard  deviation  from  the  mean. 
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August  1987  (10th  percentile,  16  observations) 
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August  1987  (10th  percentile,  64  observations) 

Latitude 


FIG.  34.  Latitudinal  variability  of  the  Tg  85  GHz  polarization  differences  at  the 
10th  percentile  threshold  for  August  1987  (groups  of  16  (top)  and  64  observations 
(bottom)).  Variability  is  ±  1  standard  deviation  from  the  mean. 
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August  1987  (60th  percentile,  all  observations) 

Latitutde 
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August  1987  (50th  percentile,  4  observations) 
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FIG.  35.  Latitudinal  variability  of  the  Tg  85  GHz  polarization  differences  at  the 
50th  percentile  threshold  for  August  1987  (all  Tb  85  observations  (top)  and  groups  of 
4  (bottom)).  Variability  is  ±  1  standard  deviation  from  the  mean. 
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August  1987  (50th  percentile,  16  observations) 

Latitude 
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August  1987  (50th  percentile,  64  observations) 
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FIG.  36.  Latitudinal  variability  of  the  Tb  85  GHz  polarization  differences  at  the 
50th  percentile  threshold  for  August  1987  (groups  of  16  (top)  and  64  observations 
(bottom)).  Variability  is  ±  1  standard  deviation  from  the  mean. 
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August  1987  (75th  percentile,  all  observations) 

Latitude 


August  1987  (76th  percentile,  4  observations) 
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FIG.  37.  Latitudinal  variability  of  the  Tg  85  GHz  polarization  differences  at  the 
75th  percentile  threshold  for  August  1987  (all  Tb  85  observations  (top)  and  groups  of 
4  (bottom)).  Variability  is  ±  1  standard  deviation  from  the  mean. 
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August  1987  (75th  percentile,  16  observations) 

Latitude 


August  1987  (75th  percentile,  64  observations) 
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FIG.  38.  Latitudinal  variability  of  the  Tb  85  GHz  polarization  differences  at  the 
75th  percentile  threshold  for  August  1987  (groups  of  16  (top)  and  64  observations 
(bottom)).  Variability  is  ±  1  standard  deviation  from  the  mean. 
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August  1987  (90th  percentile,  all  observations) 
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August  1987  (90th  percentile,  4  observations) 
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FIG.  39.  Latitudinal  variability  of  the  Tg  85  GHz  polarization  differences  at  the 
90th  percentile  threshold  for  August  1987  (all  Tb  85  observations  (top)  and  groups  of 
4  (bottom)).  Variability  is  ±  1  standard  deviation  from  the  mean. 
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August  1987  (90th  percentile,  16  observations) 

Latitude 


August  1987  (90th  percentile,  64  observations) 

Latitude 


FIG.  40.  Latitudinal  variability  of  the  Tb  85  GHz  polarization  differences  at  the 
90th  percentile  threshold  for  August  1987  (groups  of  16  (top)  and  64  observations 
(bottom)).  Variability  is  ±  1  standard  deviation  from  the  mean. 
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The  main  difference  involves  a  shift  towards  lower  ATb  values.  On  average,  the 
shift  is  less  than  5  K.  This  implies  that  resolution  does  not  severely  impact  the 
overall  characteristics  of  the  sensitivity  distributions.  Further  decreases  of  resolution 
(Fig.  34)  yield  similar  results  because  the  sensitivity  distributions  are  not 
significantly  altered  even  by  resolutions  only  1/8  the  original.  The  hemispheric 
symmetry  remains  intact  while  the  intralatitudinal  variability  remains  fairly  constant. 
The  rrCZ  regions  display  the  smallest  intralatitudinal  variability  regardless  of  the 
resolution  while  the  subtropical  regions  show  the  largest. 

The  trend  of  decreasing  ATb’s  with  decreasing  resolution  continues  for  all 
remaining  percentile  thresholds.  This  is  not  a  new  discovery;  however,  the  graphs 
show  that,  on  average,  the  polarization  difference  decreases  no  more  than  10  K 
between  the  all-observation  cases  and  the  64-observation  cases  irrespective  of  the 
percentile  threshold. 

Finally,  three  regions  that  were  selected  for  their  differing  and  diverse 
atmospheric  regimes,  were  analyzed  to  determine  the  exact  extent  resolution 
influences  sensitivity  distributions.  For  each  region,  polarization  difference  vs. 
percentile  threshold  was  plotted  for  all  four  Tb  85  group  sizes  (total,  4,  16,  and  64 
observations)(Fig.  41)  during  August  1987. 

All  three  regions  are  located  in  the  North  Atlantic.  The  first  region  is  between 
40  and  45°N  and  3  0-3  5° W.  It  is  representative  of  the  midlatitude  regions  where  large 
differences  in  cloud  coverage  and  water  vapor  can  be  found  daily  due  to  transient 
extratropical  cyclones  migrating  with  the  polar  jet  stream.  This  weather  regime 
implies  a  large  variability  in  the  polarization  differences.  The  four  midlatitude  group 
plots  have  the  least  steepest  slopes  of  the  three  selected  regions  which  is  indicative  of 
a  large  polarization  difference  variability.  The  degradation  of  the  polarization 
differences  across  the  range  of  resolutions  is  minimal;  it  is  less  than  10  K  at  the  10th 
percentile  to  approximately  12  K  at  the  90th  percentile.  Considering  the  lowest 
resolution  profile,  moderately  sensitive  algorithms  (ATb  ^  l^K)  can  be  used  80%  of 
the  time  in  this  region. 


Percentile  Percentile  Percentile 


August  1987,  Midlatitude  box  (40-45°  N,  20-25°  W) 


August  1987,  Subtropical  high  box  (25-30°N,  20-25°W) 


August  1987,  ITCZ  box  (5-10°N,  20-25°W) 


FIG.  41.  Regional  analysis  comparing  resolution  effects  on  sensitivity 
distributions  as  a  function  of  climate  regime. 
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The  second  region  is  located  is  located  between  25  and  30°N  and  20-25°W.  This 
region  is  representative  of  the  subtropical  high  belt  where  rain  and  dense  clouds  are 
minimal  due  to  large  areas  of  subsidence.  As  a  result  of  the  semi-permanent  high 
pressure  systems  there  exists  less  atmospheric  variability  in  these  regions.  The  second 
graph,  in  Fig.  41,  shows  relatively  small  ATb  variabilities  compared  to  the  midlatitude 
plots,  indicated  by  the  larger  slopes  of  all  four  group  plots.  Also,  the  plots  are 
significantly  less  opaque  evidenced  by  their  larger  ATb  values. 

Finally,  the  degradation  of  the  polarization  differences  across  the  range  of 
resolutions  is  extremely  small;  it  is  less  than  10  K  at  every  percentile.  Because  of  the 
relatively  small  atmospheric  opacity  in  Subtropical  Highs,  the  most  sensitive 
algorithms  can  be  used  nearly  all  of  the  time. 

The  third  region  is  located  between  5-10®N  and  20-25°W).  This  region  is 
representative  of  the  tropics,  specifically,  the  ITCZ,  which  is  characterized  by  strong 
upward  motion  associated  with  low  level  convergence  of  the  northeast  and  southeast 
trade  winds.  The  existence  of  ample  water  vapor,  dense  convective  clouds,  and  heavy 
rainfall,  creates  a  relatively  opaque  atmosphere  throughout  the  ITCZ.  This  large 
opacity  is  manifested  by  relatively  small  ATb  values  in  all  four  plots  of  the  third  graph 
(Fig.  41).  CLW  and  PW  variability  is  the  smallest  of  the  three  regions  even  though 
their  absolute  quantities  are  the  greatest.  The  relatively  small  variability  of  CLW  and 
PW  is  characterized  by  extremely  steep  slopes  (as  seen  in  the  third  graph). 

Finally,  the  degradation  of  polarization  differences  is  problematic  since  it  appears 
that  none  of  the  plots  are  valid  at  the  90th  percentile.  Since  atmospheric  opacity  is 
relatively  large,  robust  algorithms  must  be  used  the  majority  of  the  time.  The 
increased  resolution  capabilities  of  the  MIMR  should  alleviate  this  situation 
somewhat. 

It  is  reasonable  to  conclude  from  this  exercise  that  the  degradation  of  the  ATb 
values  resulting  from  decreased  resolutions  is  not  significant  between  frequencies  (on 
the  order  of  10  K).  Incorporating  the  6  and  10  GHz  channels  should  not  adversely 
affect  the  quality  of  the  retrievals  for  SST  and  WS. 
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CHAPTER  VI 

DISTRIBUTION  OF  CLOUD  LIQUID  WATER 

1.  Multi-frequency  retrieval  algorithm 

The  second  objective  of  this  study  was  to  understand  and  explain  the  sensitivity 
distributions  in  terms  of  the  CLW  distributions  over  selected  oceanic  areas  containing 
nonprecipitating  clouds.  To  accomplish  this,  an  operational  CLW  algorithm  from  the 
CalibrationA'’alidation  studies  (Alishouse  et  al.  1990)  was  selected. 

Their  algorithm  is  of  the  D-matrix  type;  however,  they  used  surface-based 
measurements  to  derive  new  matrix  coefficients.  The  retrieval  algorithm  (1)  is  linear 
in  brightness  temperature  and  is  configured  for  4  channels.  It  is  in  the  form  CLW  =  ao 
+  La;  *  Tei  (i  =  1,2,  ...  7),  where  the  aj’s  are  the  matrix  coefficients  and  the  Tgi’s  are 
the  brightness  temperatures  associated  with  a  specific  channel.  The  final  algorithm  is: 

CLW  (kg/m^)  =  -3.14559  +  6.0257E  -  3  ♦  19H  -  (4.8803E  -  3  *  22V)  +  (1.9595E 
-2*37V)-(3.0107E-3  *85H)  (1) 


2.  Regional  cloud  liquid  water  distributions 

Four  regions  (corresponding  to  the  climate  zones  in  Table  2)  were  selected  to 
serve  as  test  cases  for  the  algorithm’s  output:  a  high  latitude  box  extending  from  55- 
60°N  and  20-25°W,  a  midlatitude  box  extending  from  40-45°N  and  20-25°W,  a 
subtropical  high  box  extending  from  25-30°N  and  20-25°W,  and  a  tropical  box 
extending  from  10-15®N  and  20-25° W  These  regions  were  identified  because  they 
represented  a  wide  range  of  weather  regimes. 

The  pairs  of  Tb  85  observations  and  their  corresponding  lower  channel 
observations  that  were  previously  categorized  as  belonging  to  the  absorption  branch 
were  injected  into  the  algorithm.  This  was  done  to  assure  that  retrieved  values  of 
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CLW  were  associated  with  nonprecipitating  clouds  or,  at  worst,  clouds  with  very  little 
precipitation. 

The  algorithm  output,  for  a  given  region  during  August  1987,  was  plotted  on  a 
Cartesian  grid  with  cumulative  percent  on  the  vertical  axis  and  CLW  amount  (kg/m^) 
on  the  horLrontal  axis  (Figs.  42  and  43).  Together,  the  four  graphs  give  some  insight 
into  the  distribution  and  variability  of  CLW  on  a  global  scale. 

The  high  latitude  box  (Fig.  42)  indicates  that  approximately  30  %  of  the  retrieved 
CLW  values  are  negative.  Negative  CLW  values  are  not  possible.  By  definition,  in 
clear  skies,  CLW  equals  zero.  Colton  and  Poe  (1994)  identified  the  problem  of 
retrieving  negative  CLW  values  in  clear  or  nearly  clear  skies  (clear  meaning  cloud 
free)  as  inherent  in  this  algorithm. 

Two  possible  reasons  for  this  discrepancy  exist.  First,  the  algorithm  is 
overestimating  the  contribution  of  water  vapor  to  the  retrieval.  Second,  there  is  noise 
in  the  observations  on  the  order  of  .05  kg/m^  (Table  3)  which  is  within  the  range  of 
the  resulting  negative  values.  In  either  case,  it  is  reasonable  to  assign  a  value  of  zero 
to  all  negative  retrievals  and  simply  state  that  this  region  experienced  clear  skies 
approximately  30  %  of  the  time  during  August  1987. 

This  graph  also  displays  a  median  CLW  value  of  .05  kg/m^.  Since  this  median 
value  is  within  the  noise  of  the  retrievals,  the  observations  are  consistent  with  the 
skies  being  clear  more  than  50  %  of  the  time. 

Finally,  the  variability  of  CLW  is  a  function  of  the  plotted  slopes.  The  steeper 
the  slope,  the  less  the  variability.  The  slopes  of  the  high  latitude,  midlatitude,  and 
tropical  boxes  are  similar;  they  all  show  relatively  high  CLW  variability  compared  to 
the  subtropical  box  (Fig.  43).  As  mentioned  in  previous  chapters,  the  high  CLW 
variability  of  the  midlatitudes  is  directly  attributable  to  transient  extratropical 
cyclones. 

The  midlatitude  box  (Fig.  42)  shows  a  slightly  higher  median  CLW  value  of  .07 
kg/m^.  The  variation  in  the  zero  crossings  suggest  that  there  is  approximately  a  ±  .03 
kg/m^  latitudinal  shift  in  the  calibration.  However,  with  the  exception  of  the  ITCZ 
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box,  all  remaining  regions  are  consistent  with  a  relatively  high  percentage  (  up  to 
50%  in  some  cases)  of  clear  sky.  Regions  with  extremely  thin  clouds  or  CLW  values 
<  .03  kg/m^  can  be  considered  clear  sky  for  practical  purposes. 

The  subtropical  high  box  (Fig.  43)  has  a  mean  CLW  value  similar  to  the  high 
latitude  box  but  possesses  less  variability  .  Applying  the  .03  kg/m  criterion,  over 
30%  of  this  box’s  retrieval  profile  can  be  associated  with  clear  skies.  This  statistic  is 
reasonable  since  this  region  is  dominated  by  strong  subsidence  which  inhibits  cloud 
development. 

The  nCZ  box  (Fig.  43)  contains  the  highest  amounts  of  CLW  with  a  median 
value  of  .29  kg/m^.  This  box’s  relatively  high  median  value  is  a  result  of  the 
persistent  and  extensive  deep  convective  clouds  associated  with  the  tropics.  However, 
the  median  value  should  probably  be  closer  to  .17  kg/m^  owing  to  the  apparent 
positive  bias  seen  in  the  CLW  distribution.  This  bias  is  another  example  of  the 
algorithm’s  inability  to  properly  account  for  the  effect  of  water  vapor  on  atmospheric 
opacity. 

The  branch  separation  method  is  not  completely  effective  in  the  tropics.  Taking 
into  account  the  .12  kg/m^  positive  bias,  nearly  40  %  of  the  retrieved  CLW  values 
exceed  .25  kg/m^  which  is  normally  associated  with  the  cutoff  between  precipitating 
and  nonprecipitating  clouds  (Wilheit,  personal  communication).  This  CLW  value 
corresponds  to  an  atmospheric  opacity  produced  by  approximately  a  0.5  mm/hr  rain 
rate. 

This  highlights  a  shortcoming  with  the  branch  separation  method;  even  though 
observations  reside  in  the  absorption  branch,  small,  precipitation-sized  hydrometeors 
contribute  to  the  brightness  temperature  through  absorption  and  emission.  The  Tb  85 
observations  which  are  suspect  are  those  with  ATb’s  <  10  K  where  atmospheric 
opacity  is  greatest. 

Finally,  it  is  of  interest  to  quantify  (at  least  approximately)  the  magnitude  of  the 
CLW  effect  on  the  ATg  85  sensitivity  distributions.  For  example,  given  a  typical 
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midlatitude  CLW  value  of  0.07  kg/m^  (reference  Fig.  42),  we  can  quantify  the 
reduction  in  the  ATg  value  as  compared  to  a  perfectly  clear  (no  CLW  or  PW)  sky. 

The  two-component  RTM  (Fig.  7)  is  used  as  the  basis  for  the  calculations  (i.e., 
cloud  temperature  =  250  K,  0  =  53.1°,  etc.).  For  a  perfectly  clear  sky  (t  =  0),  the  ATg 
value  is  94.15  K.  After  introducing  the  0.07  kg/m  of  CLW,  and  remembering  to 
compute  for  a  two-way  pass  through  the  atmosphere,  the  optical  depth  is  .230  which 
equates  to  a  transmissivity  of  0.7942  (ATg  value  of  74.77  K).  Therefore,  CLW 
reduces  the  ATg  by  a  factor  of  approximately  1.26. 

Next,  the  all-observation  curve  for  the  midlatitude  box  in  figure  41  (which 
corresponds  to  the  midlatitude  box  in  fig.  42)  is  extrapolated  to  find  the  ATg  value  at 
the  zeroth  percentile  (38  K)  .  This  allows  us  to  compare  the  mean  value  (28  K)  of  this 
curve  against  its  clear  sky  value  (38  K).  Multiplying  the  transmissivity  with  the 
extrapolated  value  (i.e.,  0.7942  ♦  38),  yields  a  ATg  of  30.18  K  which  allows  us  to 
quantify  the  CLW’s  contribution  to  the  overall  reduction  in  polarization  difference. 

This  implies  that  the  reduction  in  the  polarization  difference  due  to  the 
contribution  of  CLW  is  approximately  93%  (28/30.18)  of  the  total.  Therefore,  in  the 
midlatitudes,  on  average,  it  is  reasonable  to  infer  that  clouds  account  for  over  90%  of 
the  variability  of  atmospheric  opacity  (at  least  during  periods  of  relatively  transparent 
atmospheric  conditions). 

The  above  procedure  was  repeated  for  the  subtropical  high  box  using  0.04  kg/m 
CLW.  The  results  indicate  that  the  CLW  contribution  to  the  reduction  in  the 
polarization  difference  is  approximately  99%  of  the  total.  Therefore,  in  the 
subtropics,  on  average,  it  is  reasonable  to  infer  that  the  clouds  account  for  nearly  all 
the  atmospheric  opacity. 

Three  of  the  four  boxes  displaying  CLW  distributions  (Figs.  42  and  43)  indicate  a 
slight  bias  towards  negative  CLW  values  on  the  order  of  .005  to  .05  kg/m  .  This 
implies  that  the  Alishouse  et  al.  algorithm  is  undercorrecting  for  the  water  vapor 
contribution.  However,  the  fourth  box  (ITCZ),  displays  a  fairly  large  positive  bias 
(approximately  .12  kg/m^)  which  implies  the  algorithm  is  overcorrecting  for  the  water 
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vapor  in  the  tropical  regions.  Therefore,  if  the  bias  (.12  kg/m  )  is  subtracted  from 
the  median  value  (.29  kg/m^),  and  the  new  median  CLW  value  of  .17  kg/m^  is  used, 
the  results  are  then  in  close  agreement  with  the  first  two  boxes.  Specifically,  clouds 
account  for  nearly  all  the  atmospheric  opacity  in  the  tropics  as  well  as  the  other 
regions. 
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CHAPTER  VII 

CONCLUSIONS  AND  RECOMMENDATIONS 

1.  Sensitivity  distribution  results 

The  global  sensitivity  distributions  reveal  that  Nelkin’s  algorithms  can  be  used 
75-90  %  of  the  time  over  the  majority  of  the  Earth’s  oceans.  This  result  varies 
depending  on  the  particular  algorithm  in  question. 

The  winter  hemisphere,  when  compared  to  the  summer  hemisphere,  provides 
better  results.  This  is  due  to  generally  dryer  and  clearer  atmospheric  conditions 
associated  with  large  scale  subsidence.  In  this  case,  the  algorithms  applicability 
reaches  90  %. 

The  sensitivity  distributions  are  affected  slightly  (approximately  10  K)  when  the 
FOV  is  increased  (resolution  decreased)  by  a  factor  of  8.  The  results  of  the  FOV 
experiment  suggest  that  the  applicability  of  the  algorithms  which  retrieve  the  ocean 
surface  parameters  SST  and  WS,  is  not  altered  significantly.  Employing  a  ATg  lower 
threshold  of  10  K  on  these  algorithms  results  in  a  decrease  in  applicability  to  a 
maximum  of  75  %  over  approximately  half  the  Earth’s  oceans. 

Conclusions  derived  from  the  algorithms  can  be  deemed  unreliable  under  various 
atmospheric  or  ocean  surface  conditions.  Certain  weather  regimes  render  the 
atmosphere  sufficiently  opaque  to  derive  unreliable  results.  These  weather  regimes 
are  characteristically  associated  with  deep  convective  clouds  with  heavy  and  persistent 
rain.  Their  locations  are  seasonally  and  latitudinally  dependent. 

A  prime  example  of  the  importance  the  time  of  year  plays  on  the  applicability  of 
the  new  algorithms  is  the  monsoon  region.  During  the  Summer,  algorithms  cannot  be 
expected  to  perform  reliably  more  than  50  %  of  the  time  whereas  during  the  Winter, 
applicability  reaches  90  %. 

The  ITCZ  and  SPCZ  are  good  examples  of  regions  where  seasonal  and  latitudinal 
factors  significantly  affect  the  global  sensitivity  distributions.  It  is  in  these  regions 
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where  robust  algorithms  will  have  to  be  selected  at  the  price  of  increased  uncertainty 
in  the  retrievals.  As  with  the  monsoon  region,  depending  on  the  location  and  intensity 
of  the  ITCZ  and  SPCZ,  which  are  seasonal  considerations,  the  algorithms  cannot  be 
expected  to  give  highly  reliable  results  more  than  50  %  of  the  time. 

Finally,  any  region  containing  sea  ice  (i.e.  Hudson  Bay  or  the  Sea  of  Okhotsk) 
cannot  be  used  by  any  of  the  4-parameter  algorithms.  Fortunately,  the  formation  of 
sea  ice  is  seasonally  dependent  and  is  limited  to  relatively  small  areas. 

2.  Cloud  liquid  water  distribution  results 

CLW  distributions  vary  with  latitude.  The  statistics  shown  in  Table  4  are  a 
compilation  of  the  regional  outputs  derived  from  the  Alishouse  et  al.  (1990)  algorithm 
for  the  North  Atlantic  during  August  1987.  Interestingly,  the  data  are  divided  into 
similar  climate  zones  as  found  in  Table  2. 

The  median  CLW  values  are  either  close  to  or  buried  within  the  noise  of  the 
observations  (.05  kg/m^).  Even  with  this  limitation,  some  knowledge  as  to  the  relative 
magnitude  of  the  amount  of  CLW  can  be  gained. 

3.  Recommendations  for  future  work 

CLW  retrievals  within  the  tropics,  especially  the  ITCZ  and  SPCZ,  contain  raining 
events.  This  is  true  even  after  the  application  of  the  branch  separation  method.  One 
possible  way  to  make  the  algorithms  more  robust  would  be  to  eliminate  Tb  85 
observations  (and  their  associated  low  channel  observations)  when  the  ATb  value  is 
less  than  10  K.  Although  this  might  lower  the  applicability  of  the  algorithms  this 
would  probably  be  offset  by  an  increased  accuracy  in  the  retrieval  amounts.  A  more 
detailed  climatology  of  the  tropics  would  need  to  be  accomplished  prior  to  obtaining 
meaningful  results. 

Additional  future  work  should  concentrate  on  the  results  obtained  from  the  global 
sensitivity  distributions.  The  idea  would  be  to  create  a  decision  tree  approach  for 
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TABLE  4.  Typical  global  cloud  liquid  water  values.  Climate  zones  based  on 
analysis  of  August  1987  regional  CLW  distributions  over  the  North  Atlantic. 


Climate  Zones 

Latitude  Range 

(Degrees) 

Median  Cloud 

Liquid  Water 
(kg/m^) 

High  Latitudes 

50-60 

.04  -  .05 

Midlatitudes 

30-50 

o 

1 

o 

Subtropics 

20-30 

.04  -  .05 

Tropics  (excluding 

0-20 

.12-  .22 

ITCZ) 
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when  to  select  a  specific  algorithm  depending  on  its  sensitivity  to  atmospheric  opacity. 
Data  could  then  be  input  into  these  algorithms,  which  are  designed  to  compute  PW 
and  other  geophysical  parameters,  and  their  output  compared  against  ground  truth  (i.e. 
radiosonde  water  vapor  profiles)  to  determine  their  validity. 
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APPENDIX  A 

SCHEMATIC  OF  THE  SSM/I’S 
CHANNELVFREQUENCY  POLARIZATIONS 


wave  propagation  (S). 

X 

Ea  also  lies  in  the  horizontal  plane  (A'). 


APPENDIX  B 

LIST  OF  ACRONYMS 


CLW 

ECMWF 

EOS 

ERT 

FOV 

ITCZ 

METOP 

MIMR 

NCAR 

PW 

RTM 

SPCZ 

SSM/I 

SST 

Tb 

ws 


-  Cloud  Liquid  Water 

-  European  Center  for  Medium-range  Weather  Forecasts 

-  Earth  Observing  System 

-  Equation  of  Radiative  Transfer 

-  Field  Of  View 

-  InterTropical  Convergence  Zone 

-  Meteorological  Operational 

-  Multi-frequency  Imaging  Microwave  Radiometer 

-  National  Center  for  Atmospheric  Research 

-  Precipitable  Water 

-  Radiative  Transfer  Model 

-  South  Pacific  Convergence  Zone 

-  Special  Sensor  Microwave/Imager 

-  Sea  Surface  Temperature 

-  Brightness  Temperature 

-  Wind  Speed 
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